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ABSTRACT 

It has been known for a long time that some red giants in globular clusters ex- 
hibit large star-to-star variations in the abundances of light elements that are not 
(N ■ 

Q ■ exhibited by field giants. This fact can be taken as evidence that the extra mix- 
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ing mechanism(s) that operate in globular cluster giants may be consequences of 
star-star interactions in the dense stellar environment. In order to constrain the 
CLi' extra mixing mechanism(s), we study the influence of helium enrichment along 

O ' the red giant branch (RGB) on the evolution of stars through the horizontal 

^ • branch. Three possible modes of helium enrichment are considered, associated 

■ with close encounters of stars in the globular clusters. We show that as a con- 
sequence of the variations in the core mass as well as in the total mass due to 

^ ' mass loss, the color of horizontal branch models are distributed over almost all 

■ range of horizontal branch. The results are discussed in relation to the scenario 
for the origin of the abundance anomalies and for the effects on the morphology 
of horizontal branch. We argue that the star-star interactions can explain not 
only the source of angular momentum of rapid rotation but also provide a mech- 
anism for the bimodal distribution of rotation rates in some globular clusters. 
We also propose the time elapsed from the latest core collapse phase during the 
gravo-thermal oscillations as the second parameter to explain the variations in 
HB morphology among the globular clusters. 

Subject headings: globular clusters: general, stars: interiors, stars: evolution, 
stars: horizontal-branch, stars: abundance, stars: rotation 
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1. Introduction 

There is growing observational evidence that points to the importance of star-star inter- 
actions in globular clusters in modifying stellar evolution, e.g., the smaller relative frequency 
of red giants in cluster cores, the overabundance of millisecond pulsars, the dependence on 
the stellar density of populations of blue stragglers (see Bailyn 1995; Hut et al. 1992, for 
reviews). In particular, the recent inflation of the number of known blue stragglers (Ferraro 
et al. 1997, 1999; Paltrinieri et al. 2001) indicates that a significant fraction of stellar pop- 
ulations may undergo encounters with neighboring stars that result in direct coalescences 
and/or binary mergers. Additional evidences may be the large star-to-star variations in the 
abundances of C, N, O, Ne, Mg and Al, which cannot be explained in terms of nucleosynthe- 
sis and mixing in the stars within the current standard framework of stellar evolution that 
takes into account only the thermal convection and molecular diffusion as the mechanisms of 
element mixing in the interior. Since these anomalous abundances are observed only among 
giants in globular clusters but not in field halo giants (see Kraft 1994; Da Costa 1997, for the 
review), it is natural to search for their origin in the environmental differences between the 
globular clusters and Galactic halo. It has already been argued that close encounters with 
nearby stars may responsible for the unusual distributions of red giants and the peculiar 
morphology of the horizontal branch at least in the central parts of globular clusters (e.g., 
Djorgovski et al. 1991; Fusi Pecci et al. 1993). 

One of characteristics of the observed abundance anomalies is that, despite large varia- 
tions in the abundances of individual elements, the sum of C, N and O abundances and the 
sum of Mg and Al abundances are found to be constant among giants in each cluster (Smith 
et al. 1996; Shetrone 1996a; Kraft et al. 1997). For C and N and also for O, the abundance 
variations are found to increase progressively along the RGB (Kraft 1994; Bellman et al. 
2001; Briley et al. 2002) to an extent which is much greater than observed for field giants 
(e.g.. Smith & Martell 2003). These facts suggest that some extra mixing mechanism or 
mechanisms operate in the course of evolution along the RGB for giants in globular clusters, 
different for field giants. The origin of the abundance anomalies is further constrained by 
the discovery of Shetrone (1996b) that, among the magnesium isotopes, ^"^Mg is depleted but 
25Mg and ^^Mg are not (see also Shetrone 1998; Yong et al. 2003). Since the burning of ^^Mg 
demands temperatures much higher than can be reached in the quiescent hydrogen burning 
shell of red giants (Langer et al. 1995), there must be machinery for raising the temperature 
in the hydrogen- containing layer high enough to burn ^^Mg. 

These observations may constitute the constraints on the nature of the mixing mecha- 
nism in cluster red giants, and from proper consideration for them, Fujimoto et al. (1999) 
propose a scenario for nucleosynthesis and extra mixing which is triggered by a close en- 
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counter with nearby stars. It has been shown by SPH simulations that during a close 
encounter, a significant amount of angular momentum can be transferred into the convective 
envelope of red giants (Davies et al. 1991, Yamada, Okazaki & Fujimoto, in preparation). 
The deposited angular momentum is transferred inwards effectively by convection to give 
rise to rotation-induced mixing in the radiative zone below the surface convective zone, and 
the latter may mix down a small amount of hydrogen from the bottom of hydrogen-burning 
shell into the upper part of helium core and causes the ignition of hydrogen shell fiashes. It 
is shown that if hydrogen reaches sufficiently deep (over a few pressure scaleheights) into the 
helium core, the hydrogen- shell flashes can be strong enough to produce the high temper- 
atures necessary to burn ^^Mg. Further, during the decay phase, matter once contained in 
the convective shell due to the hydrogen flash is dredged up by surface convection, resulting 
in the enrichment of surface layers with the nuclear products of the hydrogen shell flash. By 
computing the nucleosynthesis during the hydrogen flashes, it has been demonstrated that 
this model can reproduce the observed relationship between Na and Al abundances and the 
anti-correlation between Mg and Al abundances (Aikawa et al. 2001, 2004). 

Recently, similar abundance variations to those observed from cluster giants have been 
detected from the turn-off stars and subgiants in NGC 6752 by Gratton et al. (2001) as the 
anti-correlations between Na and O and between Mg and Al. The presence of the abundance 
anomalies in these unevolved stars has been inferred from the variations in the strengths of 
CN band and from CN and CH ant i- correlations extending from giants down to upper main 
sequence stars (Suntzeff & Smith 1991 for NGC 6752, Briley et al. 1992 for M5, Cannon et 
al. 1998 for 47 Tuc, Cohen 1999 for M71), and also from the correlation between CN band 
strength and Na abundance for 47 Tuc (Briley et al. 1996). These facts have tended to be 
interpreted as refuting the evolutionary scenario that the abundance anomalies are produced 
during the evolution along the giant branch and as favoring the primordial scenario that 
the stars were polluted by and/or born of gas already exhibiting the abundance anomalies 
(e.g., see Sneden et al. 2004). It is possible, however, to produce the abundance anomalies 
in unevolved stars through mass transfer from those giants that have already developed the 
abundance anomalies during the same kinds of star-star interactions that can deposit angular 
momentum into the convective envelope of giants (Shimada et al. 2003). Actually, the SPH 
simulations demonstrate that during a close encounter with a red giant, main sequence stars 
can accrete envelope mass from the red giant up to ~ 1O~^M0 (Davies et al. 1991, Yamada, 
Okazaki & Fujimoto, in preparation), which is enough to disguise their surface convective 
layer with the accreted material. 

For the metal-poor stars of intermediate masses, it is also shown that the temperature 
at the bottom of surface convective zone becomes high enough for ^^Mg to burn during 
the asymptotic giant branch (Ventura et al. 2001; Denissenkov & Herwig 2003). This leads 
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specifically to a recycling scenario of envelope matter ejected from erstwhile AGB stars, 
which claims that the presently observed low mass stars have been polluted by accreting 
these ejecta (Thoul et al 2002) and/or have been born from the polluted gas as second 
generation stars in the globular clusters (D'Antona 2004; D'Antona & Caloi 2004). Since 
those stars which undergo hot bottom burning also experience the third dredge-up, however, 
the pollution is expected to be necessarily attendant with enormous enhancements of CN 
abundances (Lattanzio et al. 2004; Fenner et al. 2004), and also with large enrichment of 
s-process elements, as observed in the case of CH stars. This is not the case for CN elements 
as stated above, and for the s-process elements, the contrary has been observed for the stars 
in NGC 6752 that exhibit an anti-correlation between Na and O and between Al and Mg 
(James et al. 2004). These observations impose difficulties upon the recycling scenario in 
addition to the obvious one of mass supply; since it takes more than several tenths of solar 
mass even for giants to disguise their surface with the anomalous abundances, the amount of 
mass necessary to explain the observations may well exceed the total envelope mass ejected 
from the erstwhile AGB stars that have experienced hot bottom burning. 

In contrast, the evolutionary scenario may be exempt from these difficulties. During 
hydrogen shell flashes, neither carbon nor s-process elements are produced. The surface 
disguise is responsible only for the unevolved stars, and hence, demands the order of mass in 
their surface convection (0.001 ~ O.OIM©) since the abundance anomaly can be generated 
in the stars themselves as soon as the surface convection becomes deep. The evolutionary 
scenario with star-star interactions taken into account may therefore have a better prospect 
of success and is worth further investigations. 

In the present work, we assume that deep mixing mechanism(s) operate along the RGB 
and explore their consequences on the evolution through the horizontal-branch phase. Since, 
because of a large Coulomb barrier, heavier elements such as Mg and Al can be processed only 
in the inner part of the hydrogen-burning shell, mixing responsible for Al and Mg abundance 
anomalies is likely to accompany the concomitant enrichment of surface helium (Langer & 
Hoffman 1995). The effects of possible helium mixing on the evolution have been investigated 
in relation to the morphology of the horizontal branch (Sweigart 1997a,b; Sweigart & Catelan 
1998; Caloi 2001) and to the luminosity function of the cluster giant branch (Langer et al. 
2000; Bono et al. 2001). In the existent works, helium mixing has been treated as more 
or less continuous process. This may not necessarily be true, however, since the extra 
mixing mechanism has yet to be properly established. In actuality, two different scenarios 
of helium enrichment may emerge from the models of deep mixing in globular cluster giants, 
proposed to date. One may assume that meridional currents and/or the turbulence generated 
by hydrodynamical instabilities of differential rotation carry matter from the top of the 
hydrogen-burning shell into the surface convective zone (Sweigart & Mengel 1979; Smith 



- 5 - 



& Tout 1992; Charbonnel et al. 1992; Denissenkov & Tout 2000). Or, one may assume 
that hydrogen shell flashes are involved, triggered by inward mixing of a small amount of 
hydrogen into the upper helium core presumably via the turbulence associated with the 
differential rotation, as proposed by Fujimoto et al. (1999). In the first scenario, the mixing 
may occur without interruption and can be approximated as a continuous process. In the 
second scenario, on the other hand, the mixing has to be intermittent. Although both 
types of mixing are conceivable in principle, it is only the latter mechanism that has been 
demonstrated to be able to reproduce the abundance anomalies of Al and Mg and the 
relationships between Na and Al and between Mg and Al observed in globular cluster giants 
(Aikawa et al. 2001, 2004). But it is also true that with the latter mechanism with high 
temperatures alone cannot explain the observed abundance variations of C and O much larger 
than that of Mg; there should be supplementary mechanism to mix matter processed under 
such low temperatures as C and O may burn but ^^Mg may not (Cavallo et al. 1998), which 
can be either due to the intermittent mixing with weak flashes of smaller temperatures or to 
the continuous mixing from the quiescent hydrogen shell burning, or both (Fujimoto et al. 
1999). The difference in the modes of helium enrichment associated with these extra- mixing 
mechanisms produce different effects on stellar evolution. From these latter differences we 
may be able to determine which mechanism of deep mixing nature actually prefers. 

The purpose of our work is to investigate the influences that the different modes of he- 
lium mixing may exert on the evolution of red giant and horizontal branches. In addition to 
the above-mentioned two modes of helium enrichment, we propose an additional mechanism 
that may operate after the major helium core flash is ignited at the tip of the RGB. These 
three modes are considered with particular attention to the growth in the core mass and the 
mass loss along the RGB and to the resultant variations in the luminosity and surface tem- 
perature of horizontal branch stars. In this paper, we investigate their general characteristics 
and provide a new framework for understanding the observed properties. 

In the current framework of the theory of stellar structure and evolution, one of the 
major factors that have not yet been very well incorporated is the stellar rotation. Since it 
is generally thought that the extra mixing mechanisms in the stellar interior are associated 
with rotation, the stars with abundance anomalies are also expected to provide important 
insights into the role of rotation on internal mixing (e.g., see Tassoul 2000; Pinsonneault 
1997). Indeed, for the globular clusters with large star-to-star variations of abundances such 
as M13, M92, M15, and M79, some horizontal branch stars are observed to rotate rapidly 
at velocity f rot sin i > 30 km s"^ (where i is the inclination angle of spin axis) (Peterson 
1983; Peterson et al. 1995; Cohen & McCarthy 1997; Behr et al. 2000a,b; Recio-Blanco et al. 
2002). The origin of angular momentum necessary for such rapid rotation is itself difficult to 
understand because the angular momentum ought to be effectively extracted during the main 
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sequence phase by magnetic braking (Skumanich 1972) and also during the red giant phase 
by mass loss. In fact, for turn-off stars and subgiants in some globular clusters, an upper 
bound on the rotation velocity is found as small as frot sinz < 3.5 and 4.7 km s~^ (Lucatello 
& Gratton 2003). The problem of the source of angular momentum has been addressed in 
relation to the storage of angular momentum in the stellar core by Sills & Pinsonneault (2000) 
and in relation to spin-up as a result of swallowing close planetary companion during the 
RGB evolution by Soker (1998). As criticized by Recio-Blanco et al. (2002), however, these 
scenarios suffer from defects; the first cannot specify an adequate source of large angular 
momentum in the core, and the second stumbles over the failure to detect the very existence 
of planets in an intense search in 47 Tuc by Gilliland et al. (2000). Furthermore, horizontal 
branch stars are known to display a bimodal distribution of rotation velocity, in which rapid 
rotation is restricted to the cooler horizontal branch stars [T^s < 15000°K), whereas the 
slower rotators (< 10 km s~^) spread over a wider range of effective temperature. This is 
at variance with the general expectation that faster rotation may lead to larger core and 
loss of larger envelope mass during RGB to shift the star bluer wards on HB. Consequently, 
there should be mechanism(s) not only to supply the angular momentum but also to bring 
about this bimodal distribution. In our scenario, the star-star interaction in a dense stellar 
environment is expected to be a possible source of angular momentum in rapidly rotating 
HB stars. In addition, we show that our scenario can provide an explanation also for the 
bimodal distribution of angular velocity in horizontal branch stars. 

The organization of the paper is as follows: In §2, we elaborate the modes of helium 
mixing to be considered and broadly classify them into three typical modes. In §3, we sum- 
marize the method and assumptions of the computations. The consequences of numerical 
simulations are presented in §4. Conclusions follow in §5 with a discussion, though prelimi- 
nary, on the relevance of our results to some of the observed properties of globular clusters. 
Detailed comparisons with the observations demand the specification of various parameters 
from cluster to cluster and will be discussed in a subsequent paper. 



2. Three Modes of Mixing for Surface Enrichment of Helium 

We consider three different modes of helium mixing in this paper. The first two are 
associated with the extra mixing mechanisms, suggested by other authors to work in cluster 
red giants, and differ primarily with regard to the time variations; one may regarded as 
continuous process that occurs without interrupting hydrogen shell burning and the other 
as intermittent process that recurs with hydrogen shell flashes. The last one is proposed to 
operate after the helium core flash ignites at the tip of RGB when the convective zone, driven 



- 7- 



by helium burning, extends outward into hydrogen-containing layers and engulfs hydrogen. 
We start with formulating the characteristics of three modes of helium enrichment discussed 
in this work. 



2.1. The Continuous Mixing Mode 

Continuous mixing applies to most of the extant deep mixing models which deal with 
mixing outward from the quiescent hydrogen burning shell, as studied e.g., by Sweigart & 
Mengel (1979). Mixing processes of this type may be properly described in terms of the 
strength of the mixing current and its depth of penetration into the hydrogen-burning shell 
(Boothroyd et al. 1995; Weiss et al. 2000; Denissenkov & Tout 2000). Because of the lack 
of a proper theory to constrain them, these two quantities must presently be deemed free 
parameters. Since nevertheless, the rate of helium mixing into the envelope is relevant to 
our problem, we adopt a simpler model with only one parameter that specifies the amount 
of helium produced in the hydrogen-burning shell that is carried out into the envelope. 
Denoting this parameter as /, we write the rate of change in the hydrogen abundance, X, 
in a shell as 

and the rate of change of the hydrogen and helium abundances, Xenv and Y^nv, in the envelope 

as 

dt dF ~ ~ {M - Mie)En 

Here 5h, -^h and En are, respectively, the nuclear energy-generation rate per unit mass, the 
total energy-generation rate in the star due to hydrogen burning, and the energy release per 
unit mass of hydrogen. The outer edge of hydrogen-burning shell Mie is defined as the point 
where X = X^nv — 0.01. 



It should be noted that mixing due to the meridional circulation is inhibited by a 
gradient in the mean molecular weight (Mestel 1957). From the energetics, the condition for 
mixing across layers which differ in the mean molecular weight hj Dfi is 

(Kippenhahn 1974), where fl is rotation rate at the shell, Qk (= [GM^/r^]^/^) is the local 
critical rotation rate, and Hp (= —dr/dlogP) is the pressure scale height. A similar condi- 
tion is applicable to mixing by hydrodynamical instabilities due to differential rotation, with 
the rotation rate replaced by the differential rotation rate (Fujimoto 1988). In the envelope 
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of red giants, a discontinuity in the abundance profile is formed at the first dredge-up. This 
discontinuity can be as large as AX ^ 0.05 with Hp > r/4. It has been argued, therefore, 
that the mixing of nuclear products from the hydrogen-burning shell to the surface can oc- 
cur only after the hydrogen-burning shell passes through the discontinuity unless the shell 
rotates at rates close to the local critical rotation rate and x > 0.02. This condition may 
also be applied to the mixing from the hydrogen-burning shell, across which there is a steep 
variation in mean molecular weight. Accordingly, a mixing current in the upper radiative 
zone penetrates into the hydrogen-burning shell only by a small amount, and hence, it is 
unlikely to mix an appreciable amount of helium into the envelope or to cause large varia- 
tions in heavier elements such as Mg, and Al except for very rapid rotation close to the local 
critical rate, as discussed by Sweigart & Mengel (1979). 



2.2. The Intermittent Mixing Mode 

The flash-assisted deep mixing model, proposed by Fujimoto et al. (1999), also postulates 
some extra mixing mechanism as in the other extant mixing models but it works in the 
opposite direction to the others. Instead of the direct outward transport of nuclear products 
from the hydrogen-burning shell, this model focuses on the inward transport of hydrogen 
from the bottom tail of the hydrogen-burning shell into the upper helium core and deals 
with the events consequent upon this mixing. In red giants, matter undergoes a large radial 
shrinkage while it leaves the bottom of the surface convective zone, passes through the 
hydrogen-burning shell and finally becomes incorporated in the helium core. If the matter 
has sufficient angular momentum, therefore, differential rotation ought to be enhanced to 
trigger hydrodynamical instabilities. The turbulence generated by these instabilities will 
prevail not only in layers exterior to the hydrogen burning shell but also in layers in interior 
region below the shell; these interior layers are separated from the envelope only by a very 
thin (much less than a pressure scale height) layer, and it is expected that the turbulence will 
cause material mixing between the tail of the hydrogen-burning shell and the top of helium 
core. If the turbulence is sufficiently strong, a small amount [X ~ 0.001) of hydrogen can 
be carried down into the helium core sufficiently deep (by a few pressure scale heights) in 
a sufficiently short time scale (10^ ~ 10^ yr) that eventually hydrogen shell flash will be 
ignited. Because of a small abundance of hydrogen involved and also of small ratio of the 
pressure scale height to the radial distance {Hp/r < 0.1) below the hydrogen burning shell, 
the barrier of mean molecular weight as discussed in eq. (3) is much abated in comparison 
with the outward mixing from the top of hydrogen burning shell. 

The mixing process associated with the shell flash may be outlined as follows: The 
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flash, once ignited, drives convection outward into the upper hydrogen-containing layers 
and strengthens as newly engulfed hydrogen as fuel, to extend it even beyond the site of 
the now extinct, quiescent hydrogen-burning shell. Along with the addition of entropy, the 
burning shell expands to reduce the pressure therein, and hence, the flash climbs over the 
peak stage and starts to decrease the temperature. As the hydrogen burning rate decreases, 
the flash convection retreats, and finally, the model star settles into a quiescent phase of 
stable hydrogen shell-burning again. Since the hydrogen burning shell shifts inwards during 
these events, the base of convective envelope moves inward in mass and may penetrate into 
the region that has been occupied by the flash convection to dredge up freshly synthesized 
nuclei, helium as well as the nuclear products during the quiescent hydrogen shell burning 
and during the flash under high temperatures. 

At the same time, the distribution of rotational velocity is modified both through the 
efficient transport of angular momentum by convection and through the expansion of the 
burning shell. The mixing of hydrogen is consequent upon the inward transfer of angular 
momentum via turbulent diffusion. During the following shell flash, the differential rotation 
will be largely eliminated in the shells involved in the flash-driven convective zone; the an- 
gular momentum is further carried inward, though partly carried back outward due to the 
expansion of the burning shell. After the dredge-up by the surface convection, a sharp dis- 
continuity in the rotational velocity forms at the interface between the base of the convective 
envelope and the outer edge of the erstwhile flash-driven convective zone, left behind the 
dredge-up (unless the top of the helium core and the bottom of convective envelope have 
similar angular momentum). There generated is strong differential rotation, and hence, vio- 
lent turbulence; the latter will trigger the inward mixing of hydrogen again to ignite another 
hydrogen shell flash as the hydrogen burning shell traverses the site of the former flash-driven 
convective zone and approaches to this interface. In this way, the events of flash-assisted 
deep mixing may recurs, and this epoch of recurrent flashes and mixings is expected to last 
until the initial difference in specific angular momentum between the top of the helium core 
and the bottom of the envelope convective zone is eliminated. Although the mass of the 
flash-driven convective zone may be much smaller than that of the envelope, a number of 
shell flashes during the mixing epoch may be able to produce observed variations in surface 
abundances if the initial difference in the angular momentum is sufficiently large. 

In conclusion, for the flash-assisted mixing model, the helium mixing occurs intermit- 
tently in a short time at the end of shell flashes unlike the other extant models of deep 
mixing. During the mixing epoch, the helium core may suffer from complicated changes in 
the period of alternative hydrogen flashes and quiescent burnings. At the ignition of a shell 
flash, the hydrogen-burning shell shifts inward and the core mass decreases abruptly, while 
the core mass changes very little after the ignition because of the short time-scale to burn 
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hydrogen. During the quiescent burning after the flash, the hehum core grows at the same 
rate that hydrogen burns in the shell until the next shell flash is ignited. The latter stable 
burning phase persists much longer than the flash, as estimated from the expression 

tree ~ lO^(AMfc,r/O.OOlMQ)(Xfc,r/O.l)/(LH/lO=^i:0) yr. (4) 

Here AMfc,r and X{c,r are the mass and hydrogen abundance of the hydrogen flash-driven 
convective zone, left behind the dredge-up, which may vary with the efficiency of inward 
hydrogen mixing and the strength of hydrogen shell flashes, i.e., according to how deep the 
hydrogen carried in, how far the flash-driven convection erodes the hydrogen-rich layers and 
how much mass is dredged up by the envelope convection. The resultant change in the 
core mass through the recurrent episodes in the mixing epoch depends on these factors, and 
it may happen that the core mass even decreases if the flashes are sufficiently strong so 
that the convective envelope penetrates into the site of the flash-driven convection beyond 
the erstwhile hydrogen-burning shell during the preceding stable quiescent phase. Precise 
modeling of these processes may, however, be beyond our current understandings on the 
hydrodynamical instabilities and the turbulence in the rotating stars. 

The thermal state in the core may also be affected by the behavior of the core during 
the mixing epoch. When a hydrogen shell flash ignites, the core is deprived of the upper 
shell of higher entropy and suffers from cooling at first, and, then, is heated as heat flows 
inward from the burning shell as a temperature inversion develops due to the flash. During 
the stable burning phase between flashes, the core undergoes compression at the same rate 
as hydrogen burns in the shell, which is larger as compared with standard evolution because 
of smaller hydrogen abundance in the remnant site of flash convection zone than in the 
envelope. The competition between these processes leads to a thermal balance in the core 
during the mixing epoch. It is worth noting that a similar situation occurs in thermally 
pulsating AGB wherein recurrent helium shell flashes interrupt hydrogen shell burning and 
dredge-up decreases the core mass. In the AGB case, it has been shown that the thermal 
state in the core is little affected by the intervening helium shell flashes and is determined 
primarily by the compression due to hydrogen shell burning (Fujimoto & Sugimoto 1979). 

Since the necessary equipment is an electron degenerate core and a deep convective 
envelope, flash-assisted deep mixing can work as soon as the stars evolve into red giants and 
even before the hydrogen burning shell reaches the discontinuity of hydrogen proflle produced 
by the flrst dredge-up. An additional requirement is a strong gradient in the rotation velocity 
to generate turbulence and to induce the inward mixing of hydrogen, which may in turn 
demand a large amount of angular momentum in the envelope. We may well attribute the 
latter to the deposition of angular momentum into the envelope of red giant through the 
tidal interactions at close encounters with neighboring stars in the dense environment of 
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globular clusters (Fujimoto et al. 1999). This interpretation is in accordance with the fact 
that large abundance anomalies have been observed for cluster giants but not for field giants 
in the Galactic halo (Kraft et al. 1997). This scenario also implies that the epoch of helium 
mixing may start at any stage of RGB evolution, being triggered by a close encounter. 

2.3. The Helium Flash-Driven Deep Mixing Mode 

When the helium core flash is ignited at the tip of the RGB, convection is driven 
outward by the helium burning and approaches the tail of the hydrogen-containing layer. For 
extremely metal-poor stars ([Fe/H] < —4), it is known that helium flash-driven convection 
extends through the bottom of hydrogen rich envelope and engulfs hydrogen (Fujimoto et al. 
1990, 1995, see also Fujimoto et al. 2000). For Population I and II stars, however, the outer 
edge fails to contact the hydrogen-rich envelope, but only by a few pressure scale heights, 
within the current standard framework of stellar evolution. For a M = 0.8Mq star with 
[Fe/H] = —1.5, the outer edge of the helium convective zone comes nearest, by 1.4 times a 
pressure scale- height, to the bottom of the hydrogen- rich envelope when the helium-burning 
rate decreases to Lhc = 8.8 x IO^Lq at 29.6 yr after the peak of helium burning rate of 
Lhg = 4.1 X 10^ Lq,. This shortest distance between the outer edge of the helium flash- 
driven convective zone and the tail of hydrogen containing layer is AM = 2.1 x 10"'^ Mq in 
mass, which corresponds to ~ 3 pressure scale-heights during the stationary burning phase 
just before the ignition of helium flash. 

On the other hand, we have argued in the preceding subsection for a mechanism that 
may bring about the inward mixing of hydrogen to trigger the flash-assisted deep mixing. 
Indeed, Aikawa et al. (2001) have concluded that the inward mixing of hydrogen over a few 
(2.5 ~ 3.5) pressure scale-heights is necessary to explain the abundance anomalies observed 
from M13 giants. If the same extra mixing mechanism works at the tip of the RGB, helium 
flash-driven convection can reach the hydrogen-containing layers to carry the matter inward. 
It is true that angular momentum is thought to be removed effectively from stars through 
mass loss along the RGB. And yet, the stars that experience a close encounter with another 
star at late stage of RGB evolution can retain sufficient angular momentum even at the 
tip of the RGB. If this is the case, we may expect the hydrogen injection into the helium 
flash-driven convective zone even for stars of younger populations. 

Hydrogen, thus mixed into the helium flash-driven convective zone, is further carried 
down by convection to burn and ignite a shell flash in the middle of the zone, with conse- 
quences similar to those demonstrated by HoUowell et al. (1990) for a population HI star. 
Once ignited, hydrogen shell flash develops to splits the convective zone into two zones, the 
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upper one driven by the hydrogen shell flash and the lower one by the helium shell flash. 
During the decay phase of the hydrogen shell flash, the flash convection retreats to disappear 
and as the burning shell expands, the surface convection deepens in mass to penetrate into 
the region where flash driven convection prevailed earlier. Although this sequence of events 
occurs only once, and yet, the decrease in core mass can be large enough to bring about ap- 
preciable changes in the surface abundance, as shown later. We call this dredge-up process, 
induced by hydrogen mixing into the helium-flash convection, as helium flash-driven deep 
mixing (abbreviated as He-FDDM in the following). 



3. Numerical Method and Approximations 

We compute the evolution of low mass stars along the RGB and the horizontal branch, 
taking into account the possible changes in the surface helium abundance according to the 
three modes of mixing discussed above. For the helium mixing processes, we adopt simpli- 
fying prescriptions to make them tractable by a one-dimensional stellar code. The initial 
stellar mass is set at O.SM© and the initial helium abundance ai Y = 0.24. The metallicity 
is chosen to be [Fe/H] = —1.5 with an 0.4 dex enhancement of a elements {Z = 6.8 x 10~^), 
relevant to the moderately metal-poor globular cluster M13, which has been extensively ob- 
served and is known to show large star-to-star variations in surface abundances. The initial 
CNO abundances are taken to be X12 = 7.0 x 10-^ = 2.7 x 10-^ and Xiq = 5.7 x 10"^ 
Other input physics are the same as in Iben et al. (1992) since we are interested in how the 
properties of the mixing models depend on the core mass at the time of the helium core 
flash. 

In the case of continuous mixing, we compute the variations in the abundances from 
eqs. (1) and (2) with the parameter / set constant. The mixing is started when the surface 
luminosity has reached log L/L0 = 2.1 (corresponding to core mass of Mi = 0.3058Mq) 
and is continued for the rest of the RGB life until the helium-burning rate reaches ILq. 
This choice is for the sake of simplicity and for comparison with other computations (e.g., 
Sweigart 1997a). In our model star, the discontinuity of abundance in the envelope formed 
by the flrst dredge-up is located at Mr = 0.3252Mq. 

In the case of flash-assisted deep mixing, it is of little use to follow the detailed process of 
hydrogen mixing since we currently lack the knowledge necessary to model rot at ion- induced 
mixing deflnitively. For the same reason, we ignore the change in core mass during the 
mixing epoch as well as the influence on the thermal state of the core, as discussed in the 
preceding section. Instead, we simply increase the helium abundance to a given value in a 
short timescale with the thermal state of the core little changed. Flash-assisted deep mixing 
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may start at any time on the RGB when the giant undergoes a close encounter with another 
star, as postulated in the scenario by Fujimoto et al. (1999). In order to see the dependence 
on the evolutionary stage when mixing occurs, we compute two cases with mixing occurrence 
at two different epochs: (1) the same stage [log L/Lq = 2.1) as the continuous mixing model 
starts and (2) a later stage (logL/L© = 2.6, when Mi = O.37O4M0). The time span 
between these two stages is 21.4 Myr, about two thirds of the RGB lifetime from the stage 
of luminosity logL/LQ = 2.1 to the tip of RGB in the no-mixing model. 

In the case of He-FDDM, again for the same reason, we do not follow the process of 
hydrogen engulfment. Instead, we assume simply that a sufficient amount of hydrogen is 
mixed and carried down in the helium convection at the stage of the maximum outward 
extension of helium convection to ignite hydrogen shell flash. The location of the shell 
where the hydrogen flash ignites is estimated by equating the lifetime of protons to the local 
convective turnover timescale. For obtaining the proton lifetime in the helium convective 
zone (at the stage of maximum helium convection), relevant reactions are ^^C(p, 7)^'^N and/or 
^^0(p, a)^^N, where the abundances by mass of ^^C and ^^O are set at 0.041 and 0.00068, 
respectively. Further, for simplicity, we assign a uniform distribution of hydrogen in the 
helium convective zone above the shell that hydrogen-fiash ignites, and to mimic the process 
of hydrogen engulfment, increase the abundance therein from X = 10~^ to given abundances 
of Xraix = 0.001 and 0.01 in ~ 10^ s. Then we follow the development of hydrogen shell 
flash and subsequent dredge-up by the envelope convection within the standard framework 
of stellar evolution. 



4. Results of Evolutionary Computations 

Table 1 summarizes the model parameters with the resultant surface abundance, l^nv, of 
helium, the core mass, Mi^zn, when they settle on the horizontal branch, the peak luminosity, 
Ltip, at the tip of RGB, the mass coordinate, M2,i, at the bottom of helium burning shell 
during the major helium flash, and the amount of mass loss, AM^\, through the RGB 
evolution. Here and in the following, the core is deflned as the interior to the mass shell 
where the hydrogen abundance reduces to a half of the surface abundance, i.e., X = Xcnv/2. 
These models differ only after the onset of helium mixing on the RGB. When they reach 
the base of the RGB, the surface helium abundance has already increased owing to the flrst 
dredge-up, as seen from the model without helium mixing, the difference from Fenv = 0.2491 
gives the degree of surface helium enrichment due to the deep mixing. 

In the following subsections, we flrst discuss the characteristics of RGB evolution for the 
each mixing model, and then, examine their effects on the morphology of horizontal branch 
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with the possible influence of mass loss taken into account. 

4.1. Characteristics through Red Giant Branch evolution 

4.1.1. The Continuous Mixing Models 

In the continuous mixing models, the helium mixing influences the growth rate of core 
in two opposite ways, i.e., it cut downs the helium addition to the core below the hydrogen 
burning shell while it enhances the hydrogen burning rate as a result of the decrease in 
the hydrogen abundance in the envelope. Figure 1 shows the time variations in the surface 
luminosity for various values of helium mixing rate, /. When the helium mixing starts, 
the increase of luminosity slows down to greater extent for larger rate of helium mixing /, 
which is a direct consequence of the delay of core growth by a factor of 1 — /. Then, as the 
helium abundance in the envelope augments by the helium mixing, the hydrogen burning 
rate increases, and finally, becomes larger than that of the no-mixing model at the tip of 
RGB. 

Figure 2 shows the variations in the surface helium abundance, the hydrogen burning 
rates, and the growth rates of core mass as a function of the core mass. Since we assume 
constant mixing rate /, the helium abundance, l^nv, in the envelope increases almost linearly 
with the core mass as seen from eqs. (1) and (2) ((ilog[l — Fcnv — Z] = //(I — f)d\og[M — 
Mie] with the effect of mass loss neglect). When compared at the same core masses, the 
hydrogen-burning rate is larger for larger helium mixing rate since the opacity decreases 
with the hydrogen abundance in the lower envelope adjacent to the hydrogen burning shell. 
Nevertheless, during the early stages, the growth rate of core mass remains smaller for the 
models of larger mixing rate because of larger deduction in the helium addition to the core. 
As the surface helium abundance increases, the growth rate of core rises owing both to the 
enhancement of hydrogen-burning rate and to the reduction in the nuclear energy release from 
unit envelope mass, eventually to exceed that of no-mixing model at Mi ~ 0.44 ~ 0.40Mq 
(see the analytical solution of eqs. 1 and 2 in Appendix A). Since the thermal state in the 
core is determined by the growth rate of core, the evolution to the ignition of helium core is 
consequent upon the competition between these two effects. 

Figure 3 illustrates the thermal state in the core on the density and temperature diagram 
for different mixing rates. For the models of larger mixing rates, the temperature in the 
hydrogen burning shell is higher so as to compensate the reduction in the pressure of hydrogen 
burning shell with decreasing hydrogen abundance in the envelope, which results in an 
increase in the hydrogen burning rate to meet the concomitant decrease in the opacity near 
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the burning shell. In contrast, the central temperature and entropy become lower in the 
models with larger mixing rate after the mixing is switched on, for the thermal state in the 
inner part of core ensues from the balance between the compression due to the increase in 
the core mass and the heat loss from the core. After the core mass grows more massive than 
Ml ~ O.3OM0, the temperature inversion develops because of neutrino losses from the center, 
and the maximum temperature comes to occur in the middle of core, near the shell with 
dlogn/dP \t= —1, which is defined as watershed by Fujimoto et al. (1984) since neither the 
radiation heat transfer nor the heat conduction by electrons is least efficient there. In the 
models of larger mixing rate, the maximum temperature stays smaller even after the growth 
rate of core resumes to be greater than that in the no-mixing model, as seen from the models 
with Ml = OAMq] it takes a quite while before the effects of preceding cooling in the central 
part are wiped away by the compression of core at the enhanced rate, for the inward heat 
flow from the watershed depends sensitively on the temperature and the heating of central 
part is not so efficient process. 

Consequently, the ignition of helium core flash is delayed until a larger core mass is 
reached for larger mixing rate in the range of / < 0.4, and for still larger mixing rate, the core 
mass at the helium ignition turns to decrease, as seen from the model of / = 0.6. Because of 
the cancellation between the above mentioned two effects, however, the difference in the core 
mass is rather small, as pointed out by Sweigart (1997a); in our case. AMi^zh < O.OO59M0 
and AMi^zn/AFenv < 0.0378. For larger mixing rates, the helium flash ignites in the outer 
shell because the watershed shifts outer in the mass fraction as the central density increases 
with the core mass and so does the shell of maximum temperature; the helium burning 
shell, or the bottom of helium convective zone, during the major core flash varies from 
M2i/Mi^ZH = 0.447 in the no-mixing model to 0.454 and 0.539 for the models of / = 0.20 
and / = 0.6, respectively. 



4- 1-2. The Intermittent Mixing Models 

The intermittent mixing differs from continuous mixing in that the core mass grows at 
the same rate as the hydrogen burning in the shell except for the short intermission due to 
the hydrogen shell flashes. Figures 4 and 5 show the time variations in the luminosity and 
the hydrogen-burning rate as a function of the core mass for models of intermittent mix- 
ing, respectively. The models denoted by solid and broken lines are constructed under the 
assumption that the thermal state in the core little changes during the mixing epoch with 
the latter duration artiflcially shortened, as discussed above After the mixing epoch, the 
hydrogen-burning rate, which has increased with the helium enrichment, raises the growth 
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rate of core and heats up the core above the no-mixing modeL Figure 6 illustrates the evo- 
lutionary trajectories of the central part on the density and temperature diagram. Because 
of the enhanced compression rate of core after the mixing epoch, the central temperature 
undergoes steeper rise with the central density for larger mixing rate; the whole core is kept 
hotter from the center through the hydrogen burning shell, and the helium core flash ignites 
at a smaller core and at an inner shell when compared among the models of the same mixing 
epoch. When compared among models with the same mixing rate, the heating effect of core 
is greater for an earlier occurrence of the mixing epoch, leading to the ignition of helium 
core flash at a smaller core mass. 

The tendency of the core mass at the helium ignition decreasing with the helium enrich- 
ment is in accordance with the former studies on the effects of different pristine helium abun- 
dances (e.g., Catelan et al. 1996, and references therein). The core mass Mi^zu at the helium 
ignition depends on the surface helium enrichment as AMi^zn/AFcnv = —0.19 ~ — O.ISM© 
for the models of mixing epoch at log(L/LQ) = 2.1, which may be compared with the de- 
pendence on the pristine helium abundance of dMi^zu/dY = —0.24Mq for the no-mixing 
models (Catelan et al. 1996). For a later occurrence of mixing epoch at log(L/LQ) = 2.6, the 
variation of core mass decreases to AMi zn/^^env = — O.IIM0. The models of the earlier 
mixing epoch have smaller luminosities at the tip of RGB than the model without helium 
mixing despite the surface helium enrichment. 

In the above computations, we ignore the time necessary to produce the amount of 
helium mixed into the surface for the consistency of the thermal state in the core. In order 
to assess the elongation of RGB lifetime, we compute the other extreme cases in the same 
way as in the continuous models but with / = 1 until the helium abundance in the envelope 
reaches to given values. This implies that the mixed helium is solely produced via the 
quiescent burning at a fixed core mass with concomitant increase in the helium abundance 
of envelope taken into account, and that the effects of the cooling during the mixing epoch 
is maximally taken into account. As shown in Figure 4 by dotted lines, it takes 28.3 and 6.8 
Myr to produce the amount of helium for the surface enrichment of Aycnv = 0.0793 with 
the mixing epoch at log(L/LQ) = 2.1 and 2.6, respectively. Consequently, the intermittent 
mixing models spend the RGB lifetime longer than the no-mixing model by 19.8 and 4.1 
Myr, respectively, and the elongation amounts to fairly large fraction of RGB lifetime in 
particular when the mixing epoch occurs in an early phase of RGB evolution. As compared 
with the continuous mixing model of the same helium enrichment (for which the elongation 
is 4.4 Myr), the RGB lifetime is much longer when the mixing epoch occurs earlier while 
it becomes shorter for later occurrence of mixing epoch. For the other models in Fig 4, we 
add the time necessary to produce the amount of helium mixed into the surface, estimated 
from the quiescent burning rate at a constant core mass. In these models, the reduction of 
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core mass is alleviated to be much smaller than obtained above (AMi zn/Al^nv = —0.11 
and — O.O2M0 for the mixing epoch at log L/Lq = 2.1 and 2.6, respectively). It should be 
noticed, however, that the effects of the cooing during the mixing epoch is overestimated, as 
stated in the preceding section. 

4- 1.3. The Helium Flash- Driven Deep Mixing Models 

In the present study, we postulate that an amount of hydrogen is mixed into the helium 
convective zone sufficiently to split the convection. Figure 7 shows the competition between 
the lifetime of protons and the convective turnover time scale in the helium convective zone 
as a function of mass coordinate for the stage of maximum extension of helium convection. 
Here the reaction rates are taken from Caughlan & Fowler (1988) and the convective turnover 
timescale is evaluated from the mixing length theory with the same mixing length in the 
evolution computations. The lifetime of protons decreases steeply with the depth while the 
convective turnover timescale varies slowly, and they meet to be comparable at the shell of 
Mr ~ O.3OM0 at the stage of maximum extension of helium convection. Assuming that the 
engulfed hydrogen is carried down to burn in the shell of Mr = O.SOMq, we compute two 
evolutionary sequence with artificially increasing the hydrogen abundance in the hydrogen- 
flash convection to X^ix = 0.001 and 0.01. 

Figure 8 shows the variations in the luminosity and in the hydrogen-burning rate against 
the core mass. Mi; the core is defined as interior to the shell where the hydrogen abundance 
is a half of the surface value, as stated above, and hence, the core mass varies with the 
extension of convection as well as the shell burning. When the hydrogen abundance is 
increased to X = 6.4 x 10"^ in the mixed shells, the hydrogen burning rate increases to 
Lh = 6.4 X 1O^L0 to split the convection into two zones. The upper one, driven by the 
hydrogen shell flash, develops further outward, while the lower one, driven by the helium 
shell flash, soon disappears owing to the decrease in the weight of overlying layer due to 
the hydrogen shell flash (when = 2.9 x IO^Lq). For the model of X^i^ = 0.001, the 
hydrogen shell flash reaches the peak of L^*^^ = 3.1 x IO^Lq. Even after this stage, the 
hydrogen convection continues to extend outward, and eventually reaches to the shell of 
Mr = O.49O67M0 and X = 0.748 beyond the former, now extinct, hydrogen burning shell (of 
Mr = O.49O3M0 when Lh = 9.5 x lO^L 0), as seen from this figure. The erosion of overlying 
hydrogen-rich layer raises the hydrogen abundance in the convective zone to X = 0.0016. 
Along with the subsequent expansion of hydrogen burning shell, the flash convection retreats 
finally to disappear and the hydrogen-shell burning settles in a stable quasi-equilibrium 
state at log Lu/Lq ~ 2.5. At the same time, the surface convection deepens in mass to 
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penetrates into the former site of hydrogen flash-driven convection down to Mr = 0.4745Mq 
after 2.6 x 10^ yrs from the onset of hydrogen mixing, which enriches the surface hehum 
^^env = 0.0354. This deepening of surface convection is consequent upon the inward shift 
of hydrogen-burning shell due to the hydrogen mixing and the shell flash (into the shell of 
Mr — O.SOMq). In this model, however, the hydrogen burning shell passes over the site of 
hydrogen-flash convection, left behind by the dredge-up, in relatively short time (~ 3 x 10^ 
yr) because of very small hydrogen abundance. Accordingly, the model thus enters into the 
horizontal branch with the core mass, Mi zu = 0.4745Mq. 

For the model of X^aix = 0.01, the hydrogen shell flash is slightly stronger [U^^^ = 
4.7 X IO^Lq), and drives the flash convection more outward up to the shell of O.4911M0 
(when Lff = 9.5 x 10"^Lq), as seen in Fig. 8 with the hydrogen abundance in the flash 
convection amounting to X = 0.0123. After the flash, the hydrogen shell burning settles in 
the stable quiescent state at the very bottom of the site of hydrogen flash convection since the 
nuclear timescale (rnuc = X/sh — 10^^ s) is much longer than the heat diffusion timescale 
(rdif ~ 4 X 10^° s), differently from the model of Xmix = 0.001. Accordingly, the surface 
convection reaches deeper in mass down to the shell of Mr = O.4211M0, which enriches the 
surface helium AFcnv = 0.1313. 

The distinguishing feature of this mode of mixing is the concomitant enrichment of car- 
bon and nitrogen, produced by the helium flash and processed by the hydrogen flash. Since 
the carbon abundance is Xq = 0.041 at the onset of hydrogen mixing, the CN enrichment 
in proportion to the helium enrichment as AXq^ / AY^^y ~ 0.041/(0.96 — 0.25) = 0.058. In 
our models, the carbon and nitrogen abundance in the envelope increase by factors of 26 
and 19 to Xu = 1.73 x 10"^ and Xu = 5.0 x 10"^ for the model of X^ix = 0.001, and up 
to Xi2 = 6.2 X 10~^ Xi4 = 1.845 x 10^^ for the model of X^ix = 0.01. Because of dredge- 
up, the oxygen abundance decreases slightly to Xig = 5.4 x 10""^ and Xig = 4.65 x 10""^, 
respectively. 

If the hydrogen mixing occurs at an earlier stage before the helium convection reaches its 
maximum extension, the mixed protons will burn in upper mass shells, as seen from Fig. 7; 
during the decaying phase of helium flash, both the temperature and the helium burning 
rate in the helium convective zone decline with time, and yet, since the temperature decrease 
is larger in the middle than in the bottom of convective zone, the increase in the lifetime of 
proton capture is larger than that of the convective turnover timescale. In the case of earlier 
occurrence, however, the pressure is also higher in the shell where protons burn so that the 
hydrogen flash is stronger to expel the hydrogen flash convection more outward. This entails 
larger hydrogen abundance in the hydrogen- flash convective zone, which in turn promotes 
the penetration of envelope convection when the hydrogen-burning shell settles in a stable 
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burning state. Because of this compensation, the resultant core mass will not differ so much 
when the model reaches on the horizontal branch. 



4.2. Characteristics through Horizontal Branch Evolution 

The relationship between the helium enrichment and the core mass when the model 
stars reach the horizontal branch, is summarized in Figure 9 for the three different modes 
of helium mixing. For the continuous mixing models, the core mass increases with the 
helium mixing rate for / < 0.4, and levels off or decreases for still greater mixing rate. 
The continuous mixing ends with the cores more massive than the no-mixing model and the 
largest ones among the three modes of helium mixing. In contrast, the intermittent mixing 
yields a smaller core than the no-mixing model; the core mass is a decrease function of the 
helium enrichment, and smaller for an earlier occurrence of mixing epoch. The differences 
in the core mass from the no-mixing model are much larger than for the continuous mixing 
models with the same helium enrichment. The He-FDDM gives the smallest core masses 
among the three modes of helium mixing and sets the lower limit to the core mass for a 
give helium enrichment since there is no additional production of helium to make up for the 
fraction carried out into the envelope, i.e., 

AMi.zn/AFenv > -Menv,nm/(0.96 - Fenv " AFenv) ^ "0.45 ~ -0.53 (5) 

in our cases, where Menv,nm is the envelope mass before the hydrogen mixing. 

The evolution during the central helium burning phase is exclusively determined by the 
size of helium core, and the enrichment of surface helium plays a part through the growth 
in the core mass due to hydrogen-shell burning. The mass loss through the RGB is also 
involved through the reduction of envelope mass which influences the hydrogen-burning rate 
as well as the effective temperature. In the foUowings, we first discuss the former effects, 
and then, turn into the problem of mass loss. 

4.2.1. The Luminosity and Lifetime on the HB 

Figure 10 shows the variations of the helium burning rate against the core mass interior 
to the hydrogen burning shell for the various mixing models. Since the helium burning rate 
is a monotone increasing function of the helium core mass, the continuous mixing models 
start with larger helium burning rates than the no-mixing model. For larger helium mixing 
rates, the helium burning rate is initially larger, and furthermore, increases during the core 
helium burning, to reach higher owing to the growth of core mass as well as to the reduction 
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of helium abundance in the convective core. On the other hand, the intermittent mixing 
and the He-FDDM models begin with the helium burning rates, smaller than the no-mixing 
model and decreasing for greater helium enrichment, opposite to the continuous mixing 
models. For all the modes of helium mixing, the helium core undergoes larger growth, 
and by the end of core helium burning phase, becomes more massive than for the no-mixing 
model, as a result of larger rates of hydrogen shell burning (see Figs. 2, 5, and 8). Eventually, 
the models with greater helium enrichment end with larger helium core, and in particular, 
those starting with smaller initial core masses exhibit greater increment of core mass when 
compared among those with the same helium enrichment. 

The difference in the initial helium burning rate reflects on the lifetime of horizontal 
branch phase. Figure 11 shows the time variations of mass in the hydrogen-depleted core 
through the horizontal branch for the various mixing models. Among three modes of helium 
mixing, the continuous mixing models have the shortest lives on the horizontal branch or 
the shortest duration of central helium phase because of the initially largest helium burning 
rates, which decrease rapidly with the degree of helium enrichment. On the other hand, the 
intermittent mixing and He-FDDM models may stay longer on the HB than the no-mixing 
model. The intermittent mixing first increases the HB lifetimes with the helium enrichment 
because of the decrease in the initial core mass, but it decreases for AFenv > 0.0793 since 
the enhanced hydrogen shell burning rate accelerates the growth of core to raises the helium 
burning rate and hasten the central depletion of helium. The He-FDDM models have the 
longest lifetimes on the central helium burning because of the smallest core mass at the 
beginning, which monotonically increases with the helium enrichment. Because of longer 
HB lifetime and also of the concomitant CNO enhancement, these models experience the 
largest increment of core, much larger than the intermittent mixing models with the same 
helium enrichment. 

The resultant lifetimes on the HB are summarized in Figure 12 as a function of helium 
enrichment, and compare with the lifetimes on the RG brighter than the HB (here taken 
to be a fixed value of L = 61.3Lq for the ZAHB luminosity of no-mixing model). For the 
continuous mixing, the HB lifetime is shorter than for the no-mixing model (by ~ 12% at the 
helium enrichment of AFenv = 0.0793) and decreases with the degree of helium enrichment; 
since the RGB lifetime increases with the helium enrichment, the ratio of the HB to RGB 
lifetimes decreases with the surface helium enrichment (by ~ 20% at AFcnv = 0.0793). For 
the intermittent mixing, the variation in the HB lifetime is rather small, and hence, the ratio 
between the HB and RGB lifetimes may vary greatly with the mixing epoch; for an earlier 
occurrence of log(L/L0) = 2.1, the ratio of the HB to RGB lifetimes at AFenv = 0.0793 
is smaller by ~ 33% at AFonv = 0.0793 than that for the no-mixing model, while the 
difference decreases to ~ 10% for a later occurrence of log(L/L0) = 2.6. In comparison with 
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the continuous mixing models, the intermittent mixing models give either smaller or larger 
ratios according to the difference in the mixing epoch. Lastly, the He-FDDM models have the 
largest HB lifetimes, and hence, the largest ratios among all the mixing models and the no- 
mixing model. These ratios may have relevance to the observations of R-parameter, defined 
as the number ratio between the HB stars to the RGB stars brighter than the former (Iben 
et al. 1969); for comparisons, we have to take into account the difference in the luminosity 
of horizontal branch since the duration of RGB sensitively depends on the assignment of the 
lowest luminosity, as will be discussed later. 

All the helium enriched models become brighter on the horizontal branch than the no- 
mixing model, including those with smaller cores, as seen from Figs. 1, 4, and 8, since the 
hydrogen shell-burning rate increases with the helium enrichment in the envelope and domi- 
nates the luminosity. The dependence of luminosity on the horizontal branch is A log Lhb — 
(2.1 ~ 2.2)Aycnv for the continuous mixing models, and slightly weaker A log Lhb — (1-6 ~ 
2.0)Ay^nv for the intermittent mixing models because of smaller core. The He-FDDM mod- 
els have stronger dependence A log Lhb — 2.3 and 2.7AFcnv because of the enhancement of 
CNO abundances with the helium enrichment. For proper comparisons with the observa- 
tions, however, we should take into account the change in the envelope mass due to the mass 
loss during the RGB evolution since the hydrogen burning rate depends on the envelope 
mass. 



4-2.2. The effect of mass loss 

It has been known for long time that the observed horizontal branch stars require mass 
loss, and further, that variations in the mass loss rate may be counted for the second pa- 
rameter necessary to explain the morphology of horizontal branch (Cast 1969; Iben & Rood 
1970; Rood 1973). It is true, however, that a reliable theory of mass loss is not yet avail- 
able to predict the reduction of envelope mass during the RGB phase (e.g., Catelan 2000). 
Furthermore, in applying to the stars in the globular clusters, we have to take into account 
mass loss, or mass transfer, during the close encounters with other stars, as we have pointed 
out already, which may vary from star to star and from cluster to cluster according to their 
physical conditions. Here we here start by adopting the Reimers formula (Reimers 1975) as 
a thumb rule to see an averaged effect of mass loss, i.e., 

M = 4 X 10-^^^^^^^^(M,/yr). (6) 

We compute the evolution of no-mixing model, taking account of mass loss according to 
the above formula with 77 = 1/3 (Reimers 1977), to find that the total mass reduces to 
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M = 0.6770Mq by the onset of helium flash. Apart from the reduction in the envelope 
mass, the mass loss little affects the evolution through the RGB, and in particular, leaves 
intact the core mass at the tip of RGB and the progress of helium core flash. 

Figure 13 illustrates the variations of total mass due to the mass loss against the core 
mass, computed from the luminosity and radius of the models without the mass loss; and the 
amount of the envelope mass, AM^i, lost during the whole RGB evolution is listed in the last 
column of Table 1. The mass loss becomes important as the star climbs the RGB and expands 
since the Reimers formula (6) implies the amount of mass loss in proportion to the surface 
radius if we take into account the fact that the duration of an evolutionary stage (while the 
core mass increases by a given amount) is in inverse proportion to the nuclear burning rate. 
The continuous mixing promotes the mass loss so that the integration of lost mass through 
the tip of RGB increases for larger helium enrichment, as discussed by Sweigart (1997a), 
because of larger radii reached at the tip of RGB along with the elongation of RGB lifetime. 
The intermittent mixing, however, not necessarily shows the enhancement of mass loss in 
correlation with the degree of helium enrichment. Instead, the intermittent mixing causes 
large variations in the amount of mass loss with the timing of mixing epoch. As mixing epoch 
occurs later on the RGB at larger radii, the amount of envelope mass lost during the period 
necessary to produce the helium mixed into the envelope is larger. Between our models with 
two different mixing epochs, the amount of envelope mass ejected during the mixing epoch 
increases from O.O13M0 to O.O23M0 for the helium enrichment of Al^nv = 0.0793. Further 
an earlier occurrence of mixing epoch entails the ignition of helium core flash at smaller 
core mass, and hence, smaller amount of mass loss. The both effects cooperate to make a 
difference of ~ 0.02Mq in the total envelope mass ejected during the RGB between these 
two models, which may causes significant influences on the horizontal branch evolution. 
Among the models with the same epoch, on the contrary, the above two works opposite 
direction to each other with the variation in the degree of helium enrichment. For larger 
helium enrichment, the production of helium spent for the surface enrichment takes longer 
time and increases the mass loss, while the resultant larger helium abundance in the envelope 
accelerates the evolution after the mixing epoch to ignite the helium core flash at smaller core 
masses with smaller radius and reduces the mass loss. Accordingly, the mass loss is rather 
weakly dependent on the degree of helium enrichment, or even may decrease for larger helium 
enrichment. Here we evaluate the amount of mass loss by assuming that helium is produced 
as much as carried into the envelope during the mixing epoch at a constant core mass. It is 
possible that shorter amount of helium can be produced, and hence, the core mass decreases 
during the mixing epoch if the inward mixing of hydrogen is very efficient. Even in the latter 
case, the reduction of core mass remains smaller than |AMi| < Ayenv(^ ~ Mi)/0.96 and 
will not affect so much our discussion. These characteristics, though derived by using the 
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Reimers formulae, hold to all other formulae for the mass loss as long as they give the mass 
loss rate increasing more rapidly than the luminosity. 

In conclusion, the intermittent mixing and He-FDDM models can produce the variations 
in the envelope mass among the models with the same degrees of helium enrichment as large 
as those of the continuous mixing models with the different degrees of helium enrichment. In 
Fig. 9, we compare the variations in the total mass after the mass loss during the RGB and 
in the core mass at the onset of the HB as a function of the surface helium enrichment. The 
continuous mixing models enter into the horizontal branch with the envelope mass decreasing 
with the degree of helium enrichment mainly because of the increase in the total amount 
of mass loss during the RGB evolution. On the contrary, the He-FDDM models show an 
opposite of large increase in the envelope mass with the mixing degree because of the largest 
decrease in the core mass. The intermittent mixing brings about milder decrease in the 
envelope mass with the amounts between the two extremes of the continuous mixing and the 
He-FDDM, and covers over the significant range of envelope mass by changing the timing of 
mixing epoch for the same helium enrichment. 

The above results show that the continuous mixing models have the smallest initial 
envelope masses on the horizontal branch among the models of three mixing modes with the 
same helium enrichment. But it may be the artifact that we consider only the mass loss, 
as described by Reimers-type formula. In actuality, we should take into account an extra 
mass loss during the close encounter, which is postulated in the scenario for the hydrogen 
flash-assisted deep mixing and the He-FDDM, as stated in §2. Although we are concerned 
only with the relatively distant encounters, not so destructive as to strip the envelope off 
the red giants, it is shown that such tidal interactions as deposit a significant amount of 
angular momentum into the envelope of red giants are attended also with a loss of envelope 
mass up to an order of O.OIMq, including the transfer onto the companion and/or the 
ejection from the system (Davies et al. 1991, Yamada et al. in preparation, 2005). This 
additional losses of envelope may yield a spread in the envelope masses, and yet, because of 
the statistical nature of encounters, the mutual relationship between the intermittent mixing 
and He-FDDM models will be preserved with the averaged envelope masses decreasing with 
the later occurrence of mixing epoch for the intermittent mixing and the largest envelope 
masses for the He-FDDM model. 

In order to elucidate their relations to the mixing modes, we compute the evolution 
of the model stars with the mass loss taken into account by the Reimers formula with 
rj = 1/3. Although the choice of rj is rather arbitrary, these models may give an insight 
into the effects of mass loss on the horizontal branch morphology and of variations in the 
core masses, combined with the helium mixing. Figure 14 shows the trajectories on the 



-24- 



HR diagram after the models reach the horizontal branch for the different modes of helium 
mixing, and the duration of each part of their trajectories can be read from the time variations 
of effective temperature in Figure 15. As helium is burnt in the convective core, the star first 
shrinks to raise the effective temperature owing to the increase in the central mean molecular 
weight, and then, expands to lower the effective temperature as the helium abundance in the 
convective core decreases significantly {Yc ~ 0.24 for the no-mixing model) and the central 
temperature increases (Sugimoto & Fujimoto 2000). As helium is exhausted in the center, 
the shell burning starts to expand the helium shell in the core and dwindle the hydrogen 
shell burning, contracting the star again to raise the effective temperature (Hayashi et al. 
1962), and then, the star expands as the mass in the helium-depleted core increases to evolve 
toward the asymptotic giant branch. 

The continuous mixing model starts the evolution of central helium burning from the 
bluest position than the no-mixing model and the other mixing models of the same AFcnv 
because of the smallest envelope mass; for greater helium enrichment, the initial envelope 
mass decreases monotonically (from M^m = 0.1724 to O.OIGGM© between / = 0.1 and 
0.6), and the trajectory shifts more blue-ward, as already pointed out by the earlier works 
(Sweigart 1997a; Sweigart & Catelan 1998). In addition, since the envelope mass may also 
affect the hydrogen burning rate in the shell surrounding the helium core, the continuous 
models with the mass loss reduces the luminosity on the horizontal branch significantly as 
compared with the models computed earlier without mass loss in Fig. 1. Accordingly the 
mass loss offsets the effects of increasing luminosity with the helium enrichment. For the case 
of / > 0.3 with the initial envelope mass M^nv = O.I28IM0, the luminosity tends to decrease 
despite the increase in the helium enrichment; the models of / > 0.4 with Mem = 0.0963Mq 
become even less luminous than the no-mixing model while the models without the mass 
loss are brighter by AlogL > 0.45. Smaller hydrogen burning rate also entails smaller 
growth of helium core, and longer HB lifetime, as shown in Fig. 11. As compared with the 
models without the mass loss, the growth of core is much depressed and no longer monotone 
increase function of the helium enrichment; the increment of core mass during the core helium 
burning is 5Mi = 0.0246, 0.0306, 0.0372, 0.0421, and O.O366M0 for / = 0, 0.1, 0.2, 0.3, and 
0.4, respectively, while 6 Mi = 0.0352, 0.0467, 0.0634, 0.0900, and O.133OM0 for the models 
without the mass loss. 

In contrast, the intermittent mixing models stays on redder side of the continuous models 
with the same helium enrichment. Initially, the effective temperature little depends on the 
surface helium abundance because of massive envelope (Menv ^ O.I9M0), but later, the 
models of larger helium enrichment display a larger blue-ward shift as the core mass grows to 
reduce the envelope mass. The increment in the core mass is significantly larger than those of 
continuous mixing models with the same helium enrichment because of larger initial envelope 
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mass, though much smaller as compared with the models without the mass loss {SMi = 
0.0494 and 0.0420Mo vs 6Mi = 0.0685 and 0.0643Mq for the models of different mixing 
epochs with logL/LQ = 2.1 and 2.6 and with the same helium enrichment, AFenv = 0.0793, 
computed under the assumption of the maximum cooling. See also Fig. 14). More remarkable 
is a variation in the trajectory due to the different mixing epoch, however. The model with 
the same helium enrichment reaches higher effective temperature as it has experienced the 
mixing epoch later on the RGB. 

The He-FDDM models give the lowest effective temperatures among the helium mixing 
models because of the largest reduction in the core mass due to the helium mixing; they are 
all located on redder side of the intermittent mixing models with the earlier occurrence of 
mixing epoch. For the lower helium enrichment of AFcnv = 0.0354, the model stays on the 
redder side whole through the core helium burning. On the other hand, for the larger helium 
enrichment as Al^nv = 0.1313, the model, even though starting with still smaller effective 
temperature, makes a much larger blue-ward turn in the late phase when the envelope mass 
is reduced below O.14M0; the increment of core mass is 6Mi = 0.1517Mq, and much larger 
than 6Mi = O.O5372M0 for the model of the smaller helium enrichment. 

It should be noted that the resultant envelope is subject to uncertainties in mass loss, as 
stated above. Although we find the relatively low effective temperatures of T^s < 10^ K for 
the intermittent mixing and He-FDDM models, high effective temperatures as observed can 
be readily realized if we adopt a slightly larger values for rj and/or if we take into account 
extra loss of envelope mass due to close encounters with environment stars in the dense stellar 
populations of globular clusters. The additional loss of envelope mass may cause a further 
blue-ward shift of evolutionary trajectory on HR diagram, and also, a statistical spread of 
the stars along the HB in the latter case. 

In summary, the behavior on the horizontal branch stars can differ with the mixing 
models even with the same helium enrichment because of differences in the core mass and 
in the envelope mass. As pointed out by Sweigart (1997a,b) and Langer & Hoffman (1995), 
the helium mixing causes the blue-ward shift on the horizontal branch. Both the continuous 
mixing and intermittent mixing models share this tendency. In addition, the intermittent 
mixing gives rise to variations in the location on the horizontal branch according to the timing 
of mixing epoch, shifting stars bluer-ward as they have undergone the mixing epoch later on 
the RGB evolution. Finally, the He-FDDM keeps the stars on redder side of the intermittent 
mixing models. Accordingly, even within the limited range of the surface helium enrichment, 
it is possible to generate a wide spread in the color on horizontal branch from red to blue 
with the combination of these mixing modes. In particular, from the intermittent and the 
He-FDDM, we may expect a variation in the color of horizontal branch stars not a monotone 
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function with the timing of close encounters that the giants come across. As a star undergoes 
such close encounter as leads the hydrogen flash- assisted deep mixing in later stage on the 
RGB, the descendant has to be situated on bluer side of horizontal branch, as predicted 
from the intermittent mixing model. If, on the other hand, a star has experienced the close 
encounter so late near the RGB tip to undergo the He-FDDM, it reaches the horizontal 
branch on the redder side than predicted for the intermittent mixing models even with the 
early occurrence of mixing epoch. The angular momentum injected into the envelope of red 
giants from the orbital motion at these encounters may have relevance to the fast rotation of 
horizontal branch stars and the bimodal distribution of their rotation velocities, as observed 
from the globular clusters with the abundance anomalies, as discussed below. 



5. Conclusions and discussion 

We have formulated and classified the possible mixing mechanism(s) to enrich the surface 
helium abundance in globular red giant stars into three modes, i.e., the continuous mixing, 
the intermittent mixing (or the hydrogen- flash assisted mixing), and the helium- flash driven 
deep mixing. All modes postulate some extra mixing mechanism(s) other than the thermal 
convection and molecular diffusion; the continuous mixing requires the mixing mechanism 
that can transport processed material outward from the top of hydrogen burning shell, and 
the latter two demand the mixing mechanism(s) that can carry down hydrogen from the 
bottom of hydrogen-rich layer into the top of helium core. We have computed the evolution 
of stars of mass 0.8Mq and the metallicity [Fe/H] = —1.5 with an 0.4 dex enhancement 
of a elements through the horizontal branch with the surface helium enrichment due to the 
different mixing modes taken into account and explore how these different mixing mechanisms 
influence the evolution along the red giant branch and the horizontal branch morphology. 

Our main results are summarized as follows; 
(1) The continuous and intermittent mixing modes have different, or even opposite, effects 
on the growth of core and also on the mass loss along the RGB. The continuous mixing 
delays the ignition of helium core flash until a larger core, and hence, larger luminosity, are 
reached for greater degree of helium enrichment. On the contrary, the intermittent mixing 
ends the red giant evolution with the core masses smaller for larger helium enrichment. The 
intermittent mixing may occur at any stage on the RGB and the helium core flash tends to 
ignite at smaller core mass for the occurrence of mixing epoch in an earlier stage of RGB 
evolution. 

The both modes of mixing prolong the lifetime on RGB because of the production of 
helium used for the surface enrichment. For the continuous mixing, the total amount of mass 



-27- 



loss during the RGB, as evaluated from the Reimers formulae, increases for larger helium 
enrichment. On the other hand, for the intermittent mixing, it is rather dependent on the 
timing of mixing and only weakly on the helium enrichment; the injected mass decreases for 
an earlier occurrence of mixing, the tendency of which ensues from the competition between 
the mass loss during the period to produce the extra helium and the acceleration of core 
growth after the mixing epoch. 

(2) The helium-fiash driven deep mixing (He-FDDM) occurs at the tip of RGB during 
the helium core flash and has nothing to do with the RGB evolution. For this mixing, the 
reduction in the core mass is the largest among three modes of mixing since there is no extra 
production of helium used for the surface enrichment. 

(3) When arriving at the horizontal branch, the continuous mixing model has the smallest 
envelope mass while the He-FDDM model has the largest envelope mass, and the intermittent 
mixing model falls between these two extremes, for a given degree of helium enrichment; the 
envelope mass decreases with the helium enrichment for the continuous mixing because of the 
mass loss, while increases for the He-FDDM because of the decrease in the core mass. These 
results are based only on the mass loss given by Reimers-type formulae, but in actuality, we 
have to take into account additional mass loss due to close encounters that giants undergo 
with environment stars, in particular, for the intermittent mixing and He-FDDM models. 
The latter may reduce further the envelope mass, and in addition, give rise to a statistically 
spread distribution of horizontal branch. 

(4) In general, the luminosity and the effective temperature of horizontal branch stars show 
tendency to increase with the helium enrichment. But the most critical thing in determining 
these properties is the envelope mass when the star reaches the horizontal branch, which may 
eliminate or even cancel out the effects of helium enrichment to increase the luminosity for 
such small envelope mass that gives higher effective temperatures than the red giant branch 
(T > 10'^ K). The envelope mass ensues from the competition between the growth of core 
and the mass loss, which differs among the mixing modes, and also from the loss of envelope 
mass during close encounter of red giants with nearby stars in the dense stellar environment, 
which should play an important part. 

It has been long known that the distributions of stars on the horizontal branch vary 
greatly from cluster to cluster in particular for mildly metal-poor globular clusters as studied 
here. So far, attempts have been made to explain the spread and blue-ward shift of stars 
observed along the horizontal branch in terms of the differences in their ages and of stochastic 
variations in the amount of mass loss, and recently, helium enrichment as well as merging 
of planetary systems are proposed as their cause (e.g., see Catelan 2005). The present 
results suggest that the difference in the mixing modes and also in the mixing epoch for the 
intermittent mixing can be additional factors, which may be associated with close encounters 
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with other stars that giants may experience in the dense stellar environment peculiar to the 
globular clusters. The different timing of mixing epoch gives rise to variations both in the 
amount of mass loss along the RGB and in the core mass at the ignition of helium flash, 
resulting in variation of blue-ward shifts and spread of stars along the horizontal branch, 
indifferently of the degree of helium enrichment. 

The unique feature of the timing of mixing epoch as parameter lies in the fact that it 
can generate not monotonic variation in the location of horizontal branch stars as a function 
of the time when the giants have experienced close encounter along the RGB evolution. As 
the stars come across the close encounter and undergo the mixing epoch later along the 
RGB, they tend to have smaller envelope, and hence, shift into higher surface temperatures 
on the HB. On the other hand, if the stars suffer from a close encounter very near to the 
tip of RGB, they will undergo the He-FDDM and will be situated at the redder side of the 
horizontal branch than the intermittent mixing models. This may have the relevance to the 
distribution of rotation velocity observed along the horizontal branch of globular clusters. 
Through these close encounters, the giants may once gain a significant amount of angular 
momentum to generate intense differential rotation between the envelope and the helium 
core, which can trigger the inward mixing of hydrogen and invoke the flash-assisted deep 
mixing. After undergoing the mixing event, however, they discharge the deposited angular 
momentum effectively through mass loss during the ascent of RGB, and hence, those stars 
which experience the intermittent mixing in an early stage of RGB evolution will lose much 
of the angular momentum to be slow rotators when settling on the horizontal branch. Those 
stars which undergo the close encounter and intermittent mixing later on the RGB also get 
rid of the amount of angular momentum with larger amount of envelope mass ejected through 
the mass loss, and finally settle on bluer side of horizontal branch as slow rotators. On the 
contrary, those stars which are subjected to close encounter very near to the tip of RGB and 
undergo the He-FDDM reach the redder side of HB than the intermittent mixing models, 
without losing the angular momentum acquired through the close encounter. The latter may 
correspond to very fast spinning stars found on the redder side of HB, and we may expect the 
bimodal distribution of stellar rotation along the horizontal branch. The globular clusters 
with abundance anomalies (M13, M15, M92) show the bimodal distribution of rotation 
velocity along the horizontal branch, fast rotation of v sin i ^ SO km s~^ is found only on 
the redder side of Teg = 15, 000 K, while the slow rotation of vsini < 10 km s~^ distributes 
continuously from the red to the blue (Peterson et al. 1995; Cohen & McCarthy 1997; Behr et 
al. 2000a; Recio-Blanco et al. 2002; Behr 2003a). The fastest rotation implies fi/fix ^ 0.01 at 
the tip of RGB if we assume the conservation of angular momentum, and SPH simulations of 
close encounters of giants with the main sequence stars have demonstrated that these orders 
of angular momentum can be transferred into red giants during a single encounter (Davies 
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et al. 1991, Yamada, Okazaki, Fujimoto 2005 in preparation). Consequently, our scenario 
and the present results explain only to the source of angular momentum necessary for the 
observed fast rotation but also the bimodal distribution of globular cluster horizontal ranch 
stars. This interpretation can be tested by exploring the surface abundances of fast rotating 
HB stars since their surface abundances have to be enriched in C and N as well as in helium 
such that AXcN ^ 0.058Ayenv 

The possibility that the stellar populations in the globular clusters are subject to modi- 
fications through the interactions with densely populated, environment stars, has been long 
considered and argued for some of the imprints of dynamical evolution of clusters, left on 
constituent stars (Hills & Day 1976; Djorgovski et al. 1991; Fusi Pecci et al. 1993). It is also 
true, however, that the point has been raised that the stellar encounters, as evaluated from 
the observed physical parameters of globular clusters (e.g.. Hills & Day 1976), may not be 
numerous enough to explain the observed anomalies (e.g., Djorgovski et al. 1991; Shimada 
et al. 2003). For example, recently reported is a correlation between the number of close 
binaries observed in X-ray and the stellar encounter rates in clusters (Pooley et al. 2003). 
On the other hand, it is also pointed out that for the relative frequency of blue stragglers 
to horizontal branch stars or to the red giant, no significant dependence is found on the 
cluster parameters including the expected collision rates, apart from a mild dependence on 
the central density (Piotto et al. 2004). These facts may indicate that simple comparisons 
will not be warranted between the stellar properties and the stellar interactions inferred from 
the current physical conditions. It should be noted, however, that stars with the anomalies 
have finite lifetimes, longer than the dynamical timescales of clusters, though shorter than 
the secular evolution of clusters, and hence, that the hysteresis of dynamical evolution may 
play a role in the observed properties of clusters. It is well known that the globular clusters 
experience the gravo-thermal oscillations, as first demonstrated by Sugimoto & Bettwieser 
(1983) (Bettwieser & Sugimoto 1984), and repeat the collapse and the subsequent expan- 
sion in the relaxation timescales but with wide variations in their intervals and amplitudes 
(see e.g., Makino 1996). During the collapsed phase, the timescales of stellar interactions 
decreases and stars may go through the close encounters in the central part of high stel- 
lar density. It may be true that only a small fraction of cluster stars are involved in the 
collapsed core of high density, where the specific collision rates increase significantly. But 
as already pointed out by Sugimoto (1996), it is important to take into account the mass 
segregation as well as stellar evolutions in discussing the effects of dynamical evolution of 
clusters on individual stars. Indeed, the mass segregation may increase the proportion of 
more massive stars in the central part and enhance the number of encounters during the 
core collapse involving red giants. These close encounters may result in either tidal capture 
or fiyby. Once the tidally captured binaries are formed, the giants have to experience close 
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approaches and/or exchange encounters with environment stars at much shorter timescales 
because of the increase in the cross sections due to the resonance and exchange interactions 
over the isolated stars (see e.g., Hut et aL 1992). The latter further enhances the stellar 
interactions not only to modify the stellar characteristics and evolution but also to change 
the spatial distribution through the recoil and ejection. As a result, the appearance of stellar 
components in the clusters may differ with the degrees of core collapse and of mass segre- 
gation. In addition, it critically depends on the time elapsed from the core collapse because 
owing to the finite lifetimes of RGB and HB stars (~ 0.1 Gyr), anomalous RGB and HB 
stars formed during the core collapse disappear after finishing their lives. Blue stragglers 
are also likely to form through the same kinds of encounters, and yet, may have different 
histories from them because of rather longer timescales for their formation and lives (> 1 
Gyr). Accordingly, the differences in the elapsed time from the latest core-collapse phase 
may be counted as a candidate for the second parameter, which leads to cluster-to-cluster 
variations in the abundance anomalies, observed from red giants, in the morphology of the 
horizontal branch, and also in the population of blue stragglers. 

As for as the rotation velocity of horizontal branch is concerned, it is found that field blue 
horizontal branch stars also show bimodal distributions, similar to those observed from BHB 
stars in some clusters (Peterson et al. 1983; Kinman et al. 2000; Behr 2003b; Carney et al. 
2003). This forms a sharp contrast with the fact that the abundance anomalies are observed 
only from the cluster giants but not from the field giants. Our scenario discussed above 
implies that the different source of angular momentum has to be assigned to fast rotating 
field BHB stars. In fact, the angular momentum can be supplied through the spin- up of red 
giants in binaries due to the tidal locking of stellar rotation with the orbital motion, and the 
acceleration of rotation may also be possible through the merging of sub-stellar objects (e.g., 
see Carney et al. 2003). This surely works for the field stars in the Galactic Halo, but may 
not for stars in the globular clusters since the primordial binaries of such wide separation as 
can accommodate red giants are likely to have been disrupted through the interactions with 
dense environment stars in the globular clusters before the components evolve to giants. The 
synchronization of red giants proceeds in the timescales of stellar expansion or of evolutionary 
timescale of red giants, and hence, much more slowly than the angular momentum transfer 
due to the tidal encounters. Accordingly, this may not generate strong differential rotations 
in the interior of red giant necessary to cause the inward mixing of hydrogen and produce the 
abundance anomalies as observed from cluster giants. The binaries of separations less than 
the radius at the tip of RGB (~ 0.5 au) undergo mass transfer when the primary stars evolve 
to RGB, to form blue stragglers. In actuahty it is found that the frequency of blue stragglers 
relative to the horizontal branch stars are much higher among the field stars than in the 
globular clusters (Preston & Sneden 2000). These facts again indicate that the interactions 
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with dense environment play critical roles in differentiate the evolution of stars in globular 
clusters from the filed halo stars. Proper understandings on the stellar evolution in globular 
clusters are yet to be established with influences of environments based on the dynamical 
evolution of clusters taking into consideration. 

5.1. Implications to Observations 

In this paper, we are devoted to reveal the basic properties of the different mixing 
modes and to formulate a framework to understand the abundance anomalies as observed 
from giants and the peculiar properties as observed from the horizontal branch stars at the 
same time. We introduce the close encounters of giants with environment stars as a trigger 
of the mixing and propose the elapsed time from the last core collapsed phase during the 
gravo-thermal oscillations as the second parameter for the variations in the morphology of 
horizontal branch. Based on them, we argue the possibility that the observed cluster-to- 
cluster variations can be interpreted in terms of proper physical process in globular clusters. 

The formation of abundance anomalies in red giants in line with our scenario has been 
discussed (Fujimoto et al. 1999; Aikawa et al. 2001, 2004). As for the horizontal branch stars, 
the comparisons with the observations reduce to give a scheme representing their effective 
temperature and luminosity in terms of three parameters; the total mass, the surface helium 
enrichment, and the core mass, which we have demonstrated differ according to the mixing 
modes. Although the detailed specification of the combinations of parameters should wait for 
future works, we here briefly review the implications of our scenario to observations by use of 
the present preliminary results in the following. We address the problems of the lifetime on 
HB and on RGB and the age estimation of globular clusters. In addition, we discuss about 
the luminosity functions of globular clusters and the second parameter problem. 

Figure 16 shows the ratio of the lifetime of HB model to that of RGB model brighter 
than the luminosity level of the ZAHB as a function of degree of helium enrichment. We take 
two kinds of reference luminosities of ZAHB, one fixed at the ZAHB luminosity (L = 53.9Lq) 
taken from no mixing model with mass loss, and the other set at the ZAHB luminosity taken 
from each model. These may approximate to the so-called R-ratio, defined as the number 
ratio of the HB stars to the RGB stars brighter than the luminosity level of the ZAHB 
(Iben et al. 1969). For both of the continuous and intermittent mixing models, the R-ratios 
decrease for constant reference luminosity owing to the decrease in the lifetime on HB and 
the increase of the lifetime on RGB. The reason for the increase of the RGB lifetime with 
increasing Yenv is due to the time necessary to produce the helium transported to the surface. 
For the reference luminosities, taken from ZAHB of each model, on the other hand, the R 
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ratio increases to be little dependent on the l^nv since the RGB lifetime is mainly determined 
by the stage of lowest luminosity and the lifetime ratio between the RGB and HB stars results 
from the ratio of the consumption rates of the nuclear fuels near these stages. Lastly, the 
He-FDDM models have the largest R values among all the mixing models; their R-values are 
even larger than that of no-mixing model and increase with increasing helium abundance 
in the envelope even for the fixed reference ZAMS luminosity. The R-value for the cluster 
as total may be given by the relevant combination of three modes or some of them. In 
this figure, the cluster R values fall around 1.1 for a moderate amount of enrichment in the 
surface helium abundance, which is slightly smaller than those derived from observation of 
globular clusters for [Fe/H] ^ —1.5 (see for example, Fig.l of Cassisi et al. 2003). It should 
be noticed, however, that the inclusion of convective overshooting or semi-convection will 
lengthen the HB lifetimes and expected to yield larger R values. 

The helium mixing may affect the age estimation of globular clusters. Since the mixing 
models have brighter HB luminosities than those without mixing, the AV^ method, based on 
the difference in the luminosities between the HB and Turn-off models, tends to give smaller 
age estimates when the helium mixing models are applied. From the relationship between 
the HB luminosity and the envelope helium enrichment, discussed in §4.2.1, we estimate the 
maximal changes, AMy, in the absolute magnitude due to the helium mixing, as 

d lose Z/HR 

AMv = -2.5 ^ AFenv ~ -0.18 0.4 for AFenv = 0.04 ~ 0.08. (7) 

env 

In deriving the most right-hand side member, we adopt AlogLne ~ 2.0Al^nv in the models 
without mass loss. If we take into account the reduction of HB luminosity with decreasing 
envelope mass, this imposes an upper bound to the variation in the HB luminosity, and 
hence, the scatter of HB luminosity due to the helium enrichment may fall practically in the 
range of current observational errors except for large enrichment exceeding this range. 

This difference in the HB luminosity may be transferred into an age estimate, r, of 
globular cluster, derived by the method through the luminosity dependence of r(cx 
M/Li^s oc L(i-")/" with L oc M"), giving the variation, Ar, in the age estimation as 

— = ~^^Qg ^ AFenv; ~ 0.13 ~ 0.26 for Al^nv = 0.04 ~ 0.08. (8) 

r a c/i^nv loge 

In the most right-hand side member, we assume a = 3.5 according to the standard value for 
the zero-age main sequence stars. If the helium abundances of HB stars in M13 do not exceed 
those in M3 by 0.04, as pointed out by Caloi (2001), the difference in the age estimates is 
at most ~ 13 % between models with and without the helium mixing, rather small to 
discriminate it from the current observations. Furthermore the reference level of ZAHB may 
depend on the degree and frequency of helium mixing, which introduces uncertainties in 
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adopting the fiducial ZAHB to compare with observations and affects the precise estimation 
of ages in the actual observations. 

Thanks to the high resolution of telescopes, RGB luminosity function is available from 
recent observations of globular clusters (see e.g., Cho et al. 2005) to be compared with the- 
oretical models. The duration of the stay of models at certain luminosity depends on the 
growth rate of helium core mass. For continuous mixing models, the growth rate of core 
becomes smaller than that of no-mixing model after the onset of mixing while it becomes 
larger in the late stage of RGB evolution, as discussed above. Accordingly, this mode of 
mixing may lead to the enhancement in the luminosity function in the lower part of RGB 
branch. The resultant luminosity function is dependent on the starting point of mixing and 
subject to the variations in the mixing rate. Although we assume it to begin at a fixed lumi- 
nosity, it is likely that the continuous mixing is also consequent upon the injection of angular 
momentum into the envelope of giant due to the close encounter, following the intermittent 
mixing. The intermittent mixing may have a stronger tendency toward the enhanced lumi- 
nosity function in the lower luminosity since the stagnation due to the production of extra 
helium during the mixing epoch and the acceleration of core growth after the mixing epoch. 
In this case, however, we should take into account the frequency of encounter during the 
core collapse phase, which is larger for the giants of higher luminosity and of larger radius, 
and tends to shift the luminosity function toward higher luminosity. The detail comparisons 
are awaited with these effects taken into account, which may constrain the modes of mixing 
and/or their occurrence. 

A comment is due here to the effects of star-star encounters in the globular cluster 
systems on the stellar distribution. It may be thought that stellar encounter rates are to 
be peaked in the central core of globular clusters, and so their remnants. In actuality, 
however, we should pay attention to their migrations after the star-star interactions. Such 
close encounters as transfer significant amount of angular momentum into giants are likely 
to end in the formation of wide binaries and hence, because of large increase in the cross 
sections, the binaries formed are likely to undergo additional interactions with environment 
stars. Eventually, through the resonant and/or exchange interactions, RGB stars in the 
binary can gain kinetic energy and be ejected from the core, and survive to be horizontal 
branch stars, systematically shifted toward bluer in H-R diagram. It should be noted that 
because of the reduction in the mass during the encounter and due to the mass loss, the 
giants with abundance anomalies are no longer the most massive, as used to be except for 
the neutron stars, and are possibly ejected by the exchange interactions with more massive 
turn-off stars from the binaries, originally formed by tidal captures. Since the crossing time 
scale of globular clusters (~ 10^ yr) is much shorter than their lifetimes (~ 10^ yr), such 
RGB stars as have experienced close encounters may escape from the core and be distributed 
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in the whole cluster. Therefore, strong radial concentrations toward the center may not be 
necessarily detected for the RGB stars with abundance anomalies and for the HB stars, 
as observed (Bedin et al. 2000). As stated in §1, this gives the mechanism of abundance 
anomalies, observed from dwarf and sub dwarf stars. During the encounters, the giants 
suffer from an additional mass loss of envelope, which is accreted by the counterparts of 
encounters, mostly the main sequence stars. 

So far we have discussed the characteristics of each mode of mixing and the ways to 
discriminate them. In actuality, they can be cooperative and/or complementary with each 
other. The continuous mixing can readily produce the decline of carbon abundance during 
the ascent of RGB, but not the variations in the Mg and Al abundances to such large extents 
as observed. On the other hand, the intermittent mixing can bring about the abundance 
anomalies of Mg and Al, including the reduction of ^^Mg. These two processes may cooperate 
to produce the great depletions of oxygen since the intermittent mixing itself can produce 
the extents of abundance variations similar to that of ^^Mg, which is never observed to be 
depleted to not so large extent. And yet, this does not necessarily imply that the continuous 
mixing can give rise to the surface helium enrichment; in the mildly metal-poor stars, the 
oxygen burns in the top of hydrogen burning shell in which the amount of helium produced 
remains as small as AY ~ 3Xcn — Z and it is possible to deplete oxygen with little helium 
enrichment (Sweigart & Mengel 1979). 

As for the fast rotators, it is observed that they amount to approximately one third of 
the sample of horizontal branch stars in M13 (Peterson et al. 1995). From our scenario, all 
these stars have to have experienced close encounters when they evolve near to the tip of red 
giant branch, and hence, have to occupy a fairly large fraction among the RGB stars. It is yet 
to be properly investigated how efficiently the most massive component stars assemble into 
the central part and go through the encounters, and we should wait for the future simulations 
of stellar dynamics with the mass spectra and evolution of stars appropriately taken into 
account, and hence, can follow the mass segregation and the encounters of stars including 
the tidally formed binaries. Here we point out that the collision rate with environment stars 
increases as the giants evolve to the tip, given the duration of high density due to the core 
collapse phase and also that the transfer angular momentum to the envelope is possible 
during such wide encounters as result in fly-by but not in the capture. But for sufficient 
contributions from the He-FDDM, we may seek for the other source(s) of angular momentum 
for the fast HB rotators. One of obvious candidates is the primordial binaries. In this case, 
it is necessary to reveal the conditions that such wide binaries as allow the components to 
evolve to red giants can survive to present, which varies from cluster to cluster among those 
otherwise having the similar properties. 
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A. Progress of helium enrichment for the continuous mixing 

In the continuous mixing models, the growth of helium core and the change in the 
chemical condition in the envelope due to the mixing are both related to the hydrogen 
burning rates, and hence, their rates are related as 

dMie (1 - f)LH 1-/1 C^ATenv . , . , . x . , . x 

dt X Eh Xenv f dt 

In deriving the most right-hand side member, we apply eq. (2). For constant values of /, 
this yields an integration and the change in the surface abundance is given as a function of 
core mass; given as a function of the core mass Mie] 

Xenv = 1 - l^nv " ^cnv = Xo[{M - M,,) / {M - M,,,o)Y ^ ^1 (A2) 

Here we assume that the helium mixing starts at the stage of core mass of Mie,o with a 
uniform distribution of chemical composition of X^nv = Xq and Z^nv in the envelope. The 
helium enrichment of numerical results differ slightly (< 6%) from those evaluated from this 
equation because of a gradient of helium abundance in the envelope of the initial model. 



REFERENCES 

Aikawa, M., Fujimoto, M. Y., & Kato, K. 2001, ApJ, 560, 937 
Aikawa, M., Fujimoto, M. Y., & Kato, K. 2004, ApJ, 608, 983 
Bailyn, C. D. 1995, ARA&A, 33, 133 

Bedin, L. R., Piotto, G., Zoccah, M., Stetson, P. B., Saviane, I., Cassisi, S., & Bono, G. 
2000, A&A, 363, 159 

Behr, B. B. 2003, ApJS, 149, 67 

Behr, B. B. 2003, ApJS, 149, 101 

Behr, B. B., Djorgovski, S. G., Cohen, J. G., McCarthy, J. K., Cote, P., Piotto, G., & 
Zoccah, M. 2000, ApJ, 528, 849 



- 36 - 

Behr, B. B., Cohen, J. G., & McCarthy, J. K. 2000, ApJ, 531, L37 

Bellmann, S., Briley, M. M., Smith, G. H., & Claver, C. F. 2001, PASP, 113, 326. 

Bettwieser, E., & Sugimoto, D. 1984, MNRAS, 208, 493 

Bono, G., Cassisi, S., Zoccah, M., & Piotto, G. 2001, ApJ, 546, L109 

Boothroyd, A.I., Sackmann, I.-J., & Wasserburg, G.J. 1995, ApJ, 442, L21 

Briley, M. M., Smith, V. V.,Suntzeff, N. B., Lambert, D. L., Bell, R.A., & Hesser, J. E. 1996, 
Nature,383, 604 

Briley, M. M., Smith, G. H., Bell, R. A., Oke, J. B., & Hesser, J. E. 1992, ApJ, 387, 612 
Briley, M. M., Cohen, J. G., & Stetson, P. B. 2002, ApJ, 579, L17 
Caloi, V. 2001, A&A, 366, 91 

Cannon, R. D., Croke, B. F. W., Bell, R. A., Hesser, J. E., & Stathakis, R. A. 1998, MNRAS, 
298, 601 

Carney, B.W., Latham, D.W., Stefanik, R.P., Laird, J.B., & Morse, J.A. 2003, AJ, 125, 293 
Cassisi, S., Salaris, M., & Irwin, A. W. 2003, ApJ, 588, 862 
Castellani, V., Giannone, P., Renzini, A. 1969, Ap&SS, 3, 518 
Catelan, M. 2000, ApJ, 531, 826 

Catelan, M. 2005, in APS Conf. Ser. Resolved Stellar Populations, eds. D. Valls-Gobaud & 
M. Chavez (San Francisco: APS) 

Catelan, M. de Freitas Pacheco, J.A. & Horvath, J.E. 1996, ApJ, 461, 231 

Caughlan, G. R., & Fowler, A. W. 1988, At. Data Nucl. Data Tables, 40, 283 

Cavallo, R. M., Sweigart, A. V., & Bell, R. A. 1998, ApJ, 492, 575 

Charbonnel, C, Vauclair, S., & Zahn, J.-P. 1992, A&A, 255.191 

Cho, D.-H., Lee, S.-G., Jeon, Y.-B. & Sim, K.-J. 1992, A&A, 255.191 

Cohen, J.G. & McCarthy, J.K. 1997, AJ, 113, 1353 

Cohen, J.G. 1999, AJ, 117, 2434 



-37- 



Da Costa, G. S. 1997, in Fundamental Stellar Properties: The Interaction between Obser- 
vation & Theory, lAU Symp. No. 189, eds T.R. Bedding, A.J. Booth, & J. Bavis 
(Kluwer) 193 

D'Antona, F. 2004, Mem. S. A. It., 75, 388 

D'Antona, F. & Caloi, V. 2004, ApJ,611, 871 

Davies, M. B., Benz, W., & Hills, J. G 1991, ApJ, 381, 449 

Denissenkov, P. A., & Tout, C. A. 2000, MNRAS, 316, 395 

Denissenkov, P. A., & Herwig, F. 2003, ApJ, 590, L99 

Djorgovski, S., Piotto, G., Phinney, E.S. & Chernoff, D.F. 1991, ApJ, 372, L41 

Fenner, Y., Campbell, S., Karakas, A.I., Lattanzio, J.C., & Gibson, B.K. 2004, MNRAS, 
353, 789 

Ferraro, F. R., Paltrinieri, B., Rood, R. T., & Dorman, B. 1999, ApJ, 522, 983. 

Ferraro, F. R., Paltrinieri, B., Fusi Pecci, F., Cacciari, C, Dorman, B., Rood, R. T., Buo- 
nanno, R., Corsi, C.E., Burgarella, D., & Laget, M. 1997, A&A, 324, 915 

Fujimoto, M.Y., Hanawa, T., Iben, I.Jr., & Richardson, M.B. 1984, ApJ, 278, 813 

Fujimoto, M. Y. 1988, A&A, 198, 163 

Fujimoto, M. Y. & Sugimoto, D. 1979, PAS J, 31, 1 

Fujimoto, M. Y., Iben, I. Jr. & Hollowell, D. 1990, ApJ, 349, 580 

Fujimoto, M. Y., Sugiyama, K., Iben, I. Jr. & Hollowell, D. 1995, ApJ, 444, 175 

Fujimoto, M. Y., Aikawa, M., & Kato, K. 1999, ApJ, 519, 733 

Fujimoto, M. Y., Ikeda, Y., & Iben, I, Jr. 2000, ApJ, 529, L25 

Fusi Pecci, F., Ferrano, F. R., Bellazzini, M., Djorgovski, S., Piotto, G., & Buonanno, R. 
1993, AJ, 105, 1145 

Gilliland, R. L. et al. 2000, ApJ, 545, L47 

Gratton, R. G., Bonifacio, P., Bragaglia, A., Carretta, E., Castellani, V., Centurion, M., 
Chieffi, A., Claudi, R., Clementini, G., D'Antona, F. et al. 2001, A&A, 369, 87 



- 38 - 

Hayashi, C, Hoshi, R., & Sugimoto D. 1962, Suppl. of the Prog, of Theo. Phys., 22, 1 

Hills, J.G.H & Day C.A. 1976, Astrophy. Lett. 17, 87 

Hollowell, D., Iben, I. Jr., & Fujimoto, M. Y. 1990, ApJ, 351, 245 

Hut, P., McMillan, S., Goodman, J., Mateo, M. et al. 1992, PASP, 104, 981. 

Iben, I. Jr., Rood, R. T., Strom, K. M., & Strom, S. E. 1969, Nature, 224, 1006 

Iben, I, Jr. & Rood, R.T. 1970, ApJ, 161, 587 

Iben, I. Jr. 1986, Fujimoto, M. Y., Sugimoto, D., & Miyaji, S. 1986, ApJ, 304, 217 

Iben, I, Jr., Fujimoto, M. Y., & MacDonald, J. 1992, ApJ, 388, 521 

James, G. Francois, P., Bonifacio, P., Bragaglia, A. et al. 2004, A&A, 414, 1071 

Kinman, T., Castelli, F., Cacciari, C., Bragaglia, A., Harmer, D. & Valdes, F. 2000, A&A, 
364, 102 

Kippenhahn, R. 1974, lAU Symposium 66 in Late Stages of Stellar Evolution, ed. Taylor, 
R. J. (Reidel), 20 

Kraft, R. P. 1994, PASP, 106, 553 

Kraft, R. P., Sneden, C., Smith, G. H., Shetrone, M. D., Langer, G. E., & Pilachowski, C. 
A. 1997, AJ, 113, 279 

Langer, G. E., & Hoffman, R. D. 1995, PASP, 107, 1177 

Langer, G. E., & Hoffman, R. D., Zaidins, C. S. 1997, PASP, 109, 244 

Langer, G. E., Bohe, M. & Sandquist, E. 2000, ApJ, 529, 936 

Latttanzio, J., Karakas, A., Campbell, S., Elliott, L. & Chieffi, A. 2004, Mem. S. A., It. 75, 
322 

Lucatello, S. & Gratton, R. G. 2003, A&A, 406, 691 
Makino, J. 1996, ApJ, 471, 796 
Mestel, L. 1957, ApJ, 126, 550 

Paltrinieri, B., Ferraro, F. R., Paresce, F., & de Marchi, G. 2001, AJ, 121, 3114 



- 39 - 

Peterson, R.C. 1983, ApJ, 275, 737 

Peterson, R.C. Tarbell, T.D., & Carney, B.W. 1983, ApJ, 265, 972 
Peterson, R.C, Rood, R.T. & Crocker, D.A. 1995, ApJ, 453, 214 
Pinsonneault, M. 1997, ARA&A, 35, 557 

Piotto, C, de Angeli, F. King, I. R., Djorgovski, S.G. et al. 2004, ApJ, 604,, L109 
Pooley, D., Lewin, W.H.G., Anderson, S.F., et al. 2003, ApJ, 591, L131 
Preston, G.W., & Sneden, C. 2000, AJ, 120, 1014 

Recio-Blanco, A., Piotto, C, Aparicio, A., & Renzini, A. 2002, ApJ, 572, L71 

Reimers, D. 1975, Mem. Soc. R. Sci. Liege, 6th series, 8, 369 

Reimers, D. 1977, A&A, 61, 217 

Rood, R. T. 1973, ApJ, 184, 815 

Shetrone, M. D. 1996, AJ, 112, 1517 

Shetrone, M. D. 1996, AJ, 112, 2639 

Shetrone, M. D. 1998, in Fundamental Stellar Properties, lAU Symp. No. 189, p. 158 

Shimada, M, Fujimoto, M. Y., Yamada, S, & . Sugimoto, D. 2003, in Astrophysical Super- 
computing using Particle Simulations, lAU Symp. No. 208, eds J. Makino & P. Hut 
(San Francisco: Astronomical Society of the Pacific), p. 445 

Sills, A. & Pinsonneauh, M. H., 2000, ApJ, 540, 489 

Skumanich, A. 1972, ApJ, 171, 565 

Smith, G. H. & Tout, C. A. 1992, MNRAS 256, 449 

Smith, G. H., Shetrone, M. D., Bell, R. A., Churchill, C. W., & Briley, M. M., 1996, AJ, 
112, 1511 

Smith, G. H. & Martell, S. L., 2003, AJ, 115, 1211 

Sneden, C, Kraft, R. P., Guhathakurta, P., Peterson, R. C, & Fulbright, J. P. 2004, AJ, 
127, 2162 



-40- 



Soker, N. 1998, AJ, 116, 1308 

Sugimoto, D. 1996, in Dynamical Evolution of Star Clusters — Confrantation of Theory and 
Observations, Proceedings of lAU No. 174, eds. P. Hut & J. Makino (Kluwer Academic 
Publ.), p.l 

Sugimoto, D. & Bettwieser, E. 1983, MNRAS, 204, 19 
Sugimoto, D. & Fujimoto, M.Y. 2000, ApJ, 538, 837 
Sweigart, A. V. 1997, ApJ, 474, L23 

Sweigart, A. V. 1997, in the Third Conference on Faint Blue Stars, eds. A.G.D. PhiUp, J. 
Liebert, & R.A. Saffer (Cambridge: Cambridge Univ. Press), 3 

Sweigart, A. V. & Mengel, J.G. 1979, ApJ, 229, 624. 

Sweigart, A. V. & Catelan, M. 1998, ApJ, 501, L63 

Tassoul, J-L. 2000, Stellar Rotation (Cambridge Univ. Press). 

Thoul, A., Jorissen, A., Goriely, S. Jehin, E., Magain, P., Noels, A., & Parmentier, G. 2002, 
A&A, 383, 491 

Ventura, P, D'Antona, F., Mazzitelh, I., & Gratton, R. 2001, ApJ, 550, L65 
Weiss, A. Denissenkov, P.A., & Charbonnel, C. 2000, A&A, 356, 181. 

Yong, D., Grundahl, F., Lambert, D.L., Nissen, RE. & Shtrone, M.D. 2003, A&A, 402, 985 



This preprint was prepared with the AAS IM^t;X macros v5.2. 



-41 - 



Table 1: Mixing models characteristics 



model 


/ 


Y 

^ env 


Mi,zh(M0) 


log Uk^{Lq) 


M2,(M0) 




No mixing 


0.0 


0.2491 


0.4903 


3.307 


0.2194 


0.1230 


Continuous 


0.1 


0.2850 


0.4913 


3.344 


0.2187 


0.1363 


Mixing 


0.2 


0.3284 


0.4933 


3.393 


0.2277 


0.1553 




0.3 


0.3801 


0.4945 


3.444 


0.2352 


0.1774 




0.4 


0.4438 


0.4962 


3.507 


0.2493 


0.2075 




0.6 


0.6124 


0.4960 


3.642 


0.2673 


0.2874 


Intermittent 




0.2850 


0.4833 


3.274 


0.2042 


0.1175 


Mixing 




0.3284 


0.4765 


3.293 


0.1890 


0.1153 






0.3801 


0.4706 


3.303 


0.1872 


0.1158 




1.0 


0.3284 


0.4819 


3.326 


0.2112 


0.1229 






0.3284 


0.4814 


3.323 


0.2118 


0.1315 log(L/L0)=2.6 




1.0 


0.3284 


0.4889 


3.362 


0.2417 


0.1420 log(L/L0)=2.6 


Helium Flash 


0.0 


0.2845 


0.4745 


3.307 


0.2194 


0.1230 X^i^ = 0.001 


Driven Mixing 


0.0 


0.3804 


0.4211 


3.307 


0.2194 


0.1230 Xmix = 0.01 
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Fig. 1. — Variations in the surface luminosity for the continuous mixing. Sohd hues denote 
the models with the mixing rate of / = 0.1, 0.2, 0.3, 0.4, and 0.6, from the bottom to the 
top, where the helium mixing is started at the stage of luminosity logL/L0 = 2.1. The 
dash-dotted line denotes the no-mixing model. The origin of horizontal axis is taken to be 
the age of to = 4.86 x 10^^ s from the beginning of zero-age main sequence. 
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Fig. 2. — Evolutionary characteristics as a function of core mass for the continuous mixing 
models with the mixing rate / = 0.2 (thick solid lines), 0.4(thin solid lines) and 0.6(broken 
lines) and for the no-mixing model (dash-dotted lines); the surface helium abundance, the 
hydrogen burning rate, and the growth rate of core mass are shown, respectively, in the 
panels from the top to the bottom. 
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Fig. 3. — The structure lines on the density-temperature diagram for the continuous mixing 
models with the mixing rate / = 0.2 (solid lines) and 0.6 (broken lines) and for the no-mixing 
model (/ = 0; dash-dotted lines); thick and thin lines denote the structure lines at a stage 
of core mass Mi = OAMq (thick lines) and at the stage just before the ignition of major 
helium core flash (thin lines), respectively. Two dotted lines denote those for the starting 
stages of helium mixing, log L/Lq = 2.1 and 2.6. Also plotted are the trajectories of the 
center of core up to the ignition of helium by the same lines as the structure lines. 
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Fig. 4. — Time variations in the luminosity for the intermittent mixing models; solid and 
broken curves denote models with the mixing epoch at log L/Lq = 2.1 andatlogL/L0 = 2.6, 
respectively. Dotted curves denote those which take into account the time necessary for the 
helium production carried out to the surface, and also, the maximal cooling during the 
mixing epoch. The no-mixing model is plotted by thin dash-dotted for comparison. 
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Fig. 5. — Variations in the hydrogen burning rate as a function of core mass for the inter- 
mittent mixing models; three models with the mixing epochs at log L/L© = 2.1 (solid lines) 
and a model with the mixing epochs at log L/Lq = 2.6 (broken lines). The no-mixing model 
is plotted by thin dash-dotted for comparison. 
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Fig. 6. — The same as Fig. 3 but for the intermittent mixing models with the hehum 
enrichment Y^nv = 0.3284: sohd and broken hnes denote the models with the mixing epoch 
log(L/L0) =2.1 and with log(L/L0) = 2.6, and dotted lines denotes the model with 
log(L/L0) = 2.1 and with the time to produce helium mixed into the envelope and the 
maximal cooling during the mixing epoch taken into account. The model without helium 
mixing is plotted by dash-dotted lines for comparison. 
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Fig. 7. — Lifetimes of proton, r^, in the helium convection (sohd hnes) and convective 
turnover timescales, Tconv? (broken hnes), at the stage of the maximum extension of hehum 
convection (thick hnes) and at a shghtly earher stage when the top of hehum convection 
comes within three pressure scale- heights to the hydrogen burning shell (thin lines), plotted 
against the mass coordinate M^. 
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Fig. 8. — Variations in the luminosity (thick hnes) and in the hydrogen-burning rate (thin 
hues) against the core mass for the hehum flash-driven deep mixing models with the amounts 
of mixed hydrogen, Xmix = 0.001 (solid lines) and X^i^ = 0.01 (broken lines). Also shown 
by dash-dotted lines is the no-mixing model. 
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Fig. 9. — The core mass, Mi zn (bottom thick hnes) and the total mass, M (top thin hnes) 
at the beginning of horizontal branch evolution, plotted against the helium enrichment, 
y^nv, for the three different modes of helium mixing, i.e., the continuous mixing (circle), 
the intermittent mixing with the mixing epoch at logL/L© = 2.1 (filled square) and at 
the mixing epoch at log L/Lq = 2.6 (open square) and with the maximal cooling during 
the mixing epoch taken into account (crosses), and the helium-flash driven deep mixing 
(triangles) . 
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Fig. 10. — Variations in the helium burning rate, Lhc, from the ignition of hehum core flash 
through the central and shell helium burning and to the beginning of thermal pulses on 
AGB, plotted against the mass of hydrogen-depleted core, Mi, for the continuous mixing 
models (top panel), for the intermittent mixing models (middle panel), and for the helium 
flash-driven deep mixing (bottom panel). In the middle plane, solid and broken lines denote 
the intermittent mixing models for the mixing epoch log L/Lq = 2.1 and 2.6, respectively. 
In the bottom panel, solid and broken lines indicate the models with Xmix = 0.001 and 
^mix = 0.01, respectively. Dash-dotted line in each panel denotes the no-mixing model. 
Vertical lines on the left end of each curve denote the helium core flashes and the vertical 
segments densely overlapped in the middle of each curve denote the oscillatory helium shell 
burning that is shown to occur at the beginning of helium shell burning around the helium- 
depleted core of mass M2 < O.2M0 (Iben et al. 1986). (Note that our computations do not 
take into account the overshootine of convective corel. 
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Fig. 11. — The time variations in the mass of hydrogen-depleted core for the continuous 
mixing (top panel), the intermittent (middle panel) and the helium- flash driven deep mixing 
(bottom panel). Designations of lines are the same as in Fig. 10, and symbols on each curve 
denote the end of central helium burning. For comparison, the variations in the core mass 
when the mass loss is taken into account are plotted for the continuous mixing model by 
dotted lines and for the no-mixing model by thin dash-dotted line. 
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Fig. 12. — Lifetimes on the RGB (filled symbols) and on the horizontal branch (open sym- 
bols) plotted against the degree of surface helium enrichment for three modes of helium 
mixing, i.e., the continuous mixing (circles connected with solid lines), the intermittent mix- 
ing at the mixing epoch of log(L/L0) = 2.1 (squares with broken lines) and at the mixing 
epoch of log(L/LQ) = 2.6 (larger open and filled squares), and the helium flash-driven deep 
mixing (triangles with dash-dotted lines). For the intermittent mixing models, the times 
necessary to produce the amount of helium mixed into the surface are added to the RGB 

lifetimes, estimated from the quiescent hydrogen shell burning at constant core masses (see 
T?i„ 1 1 ^ 
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Fig. 13. — The decrease in the stellar mass due to the mass loss, computed from the Reimers 
formula with 7] = 1/3, plotted as a function of mass in the hydrogen depleted core for the 
no- mixing model (dash-dotted line), for the continuous mixing models (broken lines) and 
for the intermittent mixing models with the mixing epoch of log L/Lq = 2.1 (solid lines) 
and with the mixing epoch of log L/Lq = 2.6 (dotted line), respectively. For the continuous 
mixing and the intermittent mixing with the earlier mixing epoch, two models are shown 
with different surface helium enrichment AY^nv = 0.0359 (thin lines), and 0.0793 (thick 
lines) . 
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Fig. 14. — Trajectories on HR diagram after tne models settle on the horizontal branch when 
the mass loss is taken into account according to the Reimers formula with t] = 1/3. Top panel 
for the continuous mixing models with / = 0.1, 0.2, 0,3, 0.4 from right to left: middle panel 
for the intermittent mixing models computed with / = 1.0 and AFcnv = 0.0793 at mixing 
epochs logL/Z/Q = 2.1 and 2.6 (solid and broken lines) and that of Al^env = 0.0359 (dotted 
line): and bottom panel for the helium-flash driven mixing models with X^ix = 0.001 (solid) 
and 0.01 (broken line), respectively. Dash-dotted line on each panel denotes the no-mixing 
model. 
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Fig. 16. — Ratio, R, of lifetimes on the horizontal branch to those on the red giant branch 
brighter than the reference luminosity, plotted against the degree of surface hehum enrich- 
ment for the continuous mixing (circles connected with solid lines), the intermittent mixing 
at the mixing epoch of log(L/LQ) = 2.1 (squares with broken lines), and the helium flash- 
driven deep mixing (triangles with dash-dotted lines). Open and filled symbols denote the 
values when the reference luminosity is set at the luminosities of zero-age horizontal branch 
of the no-mixing model (i.e., L = 53.9Lq) and at those of individual models, respectively; 
large open and filled squares denote those for the intermittent mixing models at the mixing 
epoch of log(L/LQ) = 2.6. The lifetimes of horizontal branch are taken from the models 
with mass loss, differently from those in Fig. 12. 



